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Section  1 
INTRODUCTION 

406 

I 

The  purpose  of  this  program  was  to  conduct  an  investigation  of  the  effects  of 

an  arbitrary  thermal  input  on  the  behavior  of  a  case -bonded  solid  propellant 

rocket  grain.  The  following  paragraphs  list  some  of  the  specific  tasks  that 

were  carried  out  under  this  program. 

(1)  Development  of  a  usable  technique  for  the  prediction  of  the  transient 
temperature  distribution  through  a  case-bonded  propellant  grain 
subjected  to  arbitrary  time -dependent  thermal  inputs.  This  was 
done  using  Biot's  approximation  solution  to  the  Fourier  heat  conduction 
equation  and  the  predictions  obtained  therefrom  were  confirmed  in  a 
series  of  motor  experiments.  The  theory  and  experimental  results 
are  presented  in  Volume  I. 

(2)  Development  of  a  usable  technique  for  the  determination  of  stress  and 
strain  in  a  case -bonded  real  viscoelastic  propellant  grain  subjected  to 
the  above  arbitrary  thermal  inputs  .  This  task  (reported  in  Volume  II) 
was  accomplished  with  Purdue  University  under  subcontract;  Purdue 
conducted  the  theoretical  study  and  Lockheed  Propulsion  Company 
carried  out  the  computer  programing  and  solution  (Volume  I)  with  the 
assistance  of  the  Mathematical  Section  of  the  Lockheed-California 
Company  and  their  IBM  7090  facility. 

(3)  Determination  of  the  pertinent  viscoelastic  and  thermal  characteriza¬ 
tions  of  propellant  materials  for  use  as  described  in  items  (l)  and  (2) 
above.  This  determination  was  made  at  LPC  for  typical  rubber -base 
propellants.  The  methods  used  and  the  results  obtained  are  reported 
in  Volume  I. 

(4)  Determination  of  applicable  failure  criteria  for  thermal  cycling  and, 
using  analog  motors,  the  confirmation  of  these  criteria.  The  results 
are  contained  in  Volume  I. 

(5)  Reduction  to  engineering  practice  of  the  above  analytical  methods,  using 
limit  analysis  as  well  as  engineering  methods  for  gravitational  slump 
and  acceleration  set-back.  These  results  are  presented  in  the  form  of 
an  engineering  handbook  (Volume  III)  which  also  includes  simplified 
methods  for  determination  of  the  viscoelastic  operators. 

(6)  Application  of  the  above  methods  for  the  determination  of  the  thermal 
environment  capabilities  of  various  large  booster  motors.  These 
results  are  contained  in  Volume  III-A,  which  is  a  classified  addendum 
to  Volume  III. 
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(7 )  Determination  of  the  existence  of  a  generalized  cyclic  thermal  input 
over  a  fixed  temperature  range  which  will  produce  a  greater  tendency 
to  failure  than  a  step  function  input  over  the  same  range.  This  problem 
was  not  resolved  analytically.  However,  a  series  of  calculations  for  a 
wide  range  of  motor  sizes  and  input  variations  was  carried  out.  Within 
the  limits  of  accuracy  of  the  computer  solution,  no  situation  was  found 
wherein  the  cyclic  input  caused  stresses  or  strains  exceeding  those 
produced  by  a  step  input.  The  worst  failure  condition  still  was  found 

to  be  that  of  a  motor  subjected  to  rapid  cooling  to  a  low  temperature, 
followed  --  after  thermal  equilibrium  was  established  --by  rapid 
external  heating.  In  one  sense  this  is  a  cyclic  input,  but  not  in  the 
sense  of  the  originally  posed  question  regarding  stress  accumulation 
beyond  the  step  input  condition.  Questions  of  cyclic  fatigue  are  not 
considered  the  same  as  stress  build-up. 

(8)  Determination  of  a  method  of  stress  analysis  using  dynamic  birefring¬ 
ence  on  viscoelasticity  --or  photoviscoelasticity,  as  it  is  frequently 
called.  The  experimental  work  was  carried  out  by  New  York  Univer¬ 
sity  and  fully  reported  in  previously  published  quarterly  reports  to  this 
program.  Unfortunately,  the  plasticized  Hysol  material  csed  was  tested 
over  a  temperature  range  which  included  the  transition  point  temperature 
of  the  material.  Further,  the  experimental  results  did  not  superpose, 
indicating  that: 

(a)  the  material  was  not  linearly  viscoelastic  or  that 

(b)  the  experimental  determinations  of  creep  and  temperature 
were  not  carried  out  with  sufficient  precision. 

Basically,  however,  it  was  determined  during  the  course  of  this  pro¬ 
gram  that  a  proper  theory  for  the  interpretation  of  stresses  using 
photoviscoelasticity  did  not  factually  exist.  This  problem  was  then 
defined  and  a  theoretical  treatment  is  given  in  Appendix  B  of  Volume  I. 
Experimental  confirmation  is  now  in  progress  under  other  urograms 
at  the  California  Institute  of  Technology,  the  University  oi  Washington, 
and  at  Lockheed  Propulsion  Company. 

Thus,  to  accomplish  the  subject  purpose  of  the  program,  and  to  present  the 
results  in  a  readily  usable  form,  the  main  body  of  this  report  has  been  divided 
into  five  main  sections. 

Section  2,  Thermal  Field  Determination,  presents  all  of  the  necessary  analyti¬ 
cal  and  experimental  techniques  for  the  determination  of  the  time -dependent 
temperature  distribution  in  a  circular  port  grain  subjected  to  an  arbitrary 
external  temperature  variation.  The  results  predicted  by  theory  are  shown 
to  be  accurately  confirmed  in  experimental  grains . 
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One  aspect  of  the  above  analysis  which  has  not  been  treated  is  that  concerning 
the  determination  of  the  heat  transfer  boundary  conductance.  Effectively,  the 
boundary  conductance  has  been  considered  herein  as  infinite;  the  motor  case 
thermal  conductivity  also  was  considered  as  infinite.  For  the  purpose,  these 
assumptions  give  answers  on  the  safe  side.  In  actual  fact,  these  assumptions 
affect  the  answer  but  slightly  for  an  insulating  or  semi -insulating  material 
such  as  rocket  grains,  and  experimental  results  confirm  this  fact.  Thus, 
using  the  methods  delineated  herein,  tne  problems  of  high  speed  aerodynamic 
induced  heat  transfer  are  readily  computed. 

The  computational  demands  upon  the  user  are  reduced  by  the  inclusion  of  a 
series  of  dimensionless  graphs  which  enable  ready  scaling  of  results  for 
various  sizes  of  rocket  motors. 

Section  3,  Propellant  Viscoelastic  Characterization,  presents  a  brief  review 
of  reduced  variable  relationships  for  viscoelastic  materials.  In  addition,  it 
presents  that  the  utilization  of  approximate  power  law  relationships  leads  to 
greatly  reduced  effort  in  the  utilization  of  the  previously  mentioned  scaling 
laws . 

Section  4,  Grain  Thermal  Stress  Analysis,  presents  the  actual  purpose  of 
this  work.  The  summary  of  results  is  presented  in  this  volume  and  the  com¬ 
plete  details  of  the  analysis  are  given  in  Volume  II.  The  simplified  limit 
analysis  is  given  whereby  the  maximum  stress  or  strain  for  a  given  tempera¬ 
ture  difference  may  be  calculated.  For  temperature  cycling,  the  full  solution 
is  required.  Some  of  the  examples  show  that  the  maximum  stress  unaer  a 
cycled  condition  may  be  considerably  less  than  the  limit  solution  --a  result 
that  may  be  extremely  useful  for  marginal  grain  designs  .  Section  4  includes 
a  discussion  of  the  theoretical  aspects  of  photothermoviscoelasticity. 

Section  5,  Grain  Failure  Analysis,  concerns  the  experimental  verification  of 
the  pertinent  failure  criteria  associated  with  the  thermal  problems  treated. 

The  two  appendices  to  this  volume  treat  more  fully  the  subjects  of  Thermal 
Field  Solutions  (Appendix  A),  and  the  Theory  of  Photovis  c  o  e  1  a  s  t  i  c  i  t  y 
(Appendix  B). 
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Section  2 

THERMAL  FIELD  DETERMINATION 


2.1  GENERAL 

Thermal  stress  determination  in  elastic  solids  requires  knowledge  only  of  the 
instantaneous  spatial  distribution  of  temperature.  Furthermore,  there  are 
many  thermo -elastic  problems  requiring  only  knowledge  of  the  average  tem¬ 
perature  for  approximate  solutions  to  the  stress -strain  problem. 

As  might  be  expected,  however,  viscoelastic  solids  require  not  only  the 
instantaneous  spatial  temperature  field  determination,  they  also  require  the 
determination  of  the  entire  time  history  of  the  spatial  temperature  field. 

Thus,  this  section  is  devoted  to  the  following; 

•  Analytical  determination  of  time  varying  spatial  temperature 
distributions  (Section  2.2.1  and  2.2.2) 

•  Numerical  solution  of  the  resulting  thermal  field  problem 
for  arbitrary  temperature  inputs  (Section  2.2.2) 

•  Experimental  determination  of  the  pertinent  propellant 
and  motor  thermal  properties  (Section  2.3) 

•  Experimental  verification  of  the  temperature  distribution 
with  the  above  thermal  properties  in  cylindrical  motor 
grains,  and  in  circular  port  motors.  (Section  2.4) 

The  laboratory  determination  of  the  pertinent  propellant  thermal  properties 
(density,  p;  conductivity,  c  ;  and  specific  heat, k)  indicated  a  rather  trivial 
variation  of  these  properties  with  temperature  over  the  temperature  range 
of  interest,  -100°F  to  +  IbO^F.  Thus,  to  make  the  problem  more  tractable  to 
analysis,  the  assumption  is  made  throughout  that  the  above  thermal  properties 
do  not  vary  with  temperature. 


2-  i 


LOCKHEED  PNOHULSION  COMPANY 


578-F 


Again,  with  the  intent  of  producing  a  tractable  analytical  formulation  of  the 
problem,  Biot's  approximate  solution  to  the  Fourier  heat  conduction  equation 
is  used  rather  than  attempting  a  rigorous  solution.  (Boley  and  Weiner,  I960, 
p-224f) 

The  excellent  correlation  between  the  predicted  temperature  fields  in  both 
time  and  space  and  the  resulting  measurement  in  test  motors  justifies  the  use 
of  all  the  aforementioned  simplifying  assumptions. 
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2.2  TEMPERATURE  FIELDS  IN  CIRCULAR  PORT  GRAINS 
2.2.  1  Biot's  Approximation  Method 


The  solution  of  the  Fourier  heat  conduction  equation  usually  leads  to  forms  of 
solutions  which  are  computationally  cumbersome.  As  an  engineering  alter¬ 
native,  therefore,  Biot's  approximation  method  (Boley  and  Weiner,  I960, 
p.226)  will  be  used.  A  brief  description  is  given  below  for  the  case  of  axi- 
symmetrical  temperature  input  to  the  outside  of  an  infinitely  long  cylindrical 
propellant  grain  with  temperature  independent  thermal  properties  and  an  in¬ 
sulated  circular  port. 

The  temperature,  T,  must  satisfy  the  heat  conduction  equation  (one -dimensional 
because  of  axial  symmetry) 


k 


pc 


dT 

TT 


where 

T 

n 

t 

k 

P 

c 


=  temperature 
=  direction  of  heat  flux 
=  time 

=  specific  heat 
=  density 

=  thermal  conductivity 


and  defining  a  heat  flow  vector  H  by 
d^T  d  T  d^  H 

so  that  it  follows 
d  H 

pc  T  =  -  ~  or  H  =  -  J  pcTdn 

dn  q, 

and 

dT  _ 

^  dn  '  ■  dt 


where  qj  is  the  depth  to  which  heat  has  penetrated  at  time,  t. 


(2-1) 


(2-2) 


(2-3) 


(2 -3a) 
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Employing  variational  techniques,  one  can  derive  functions 
V  =  Y  p  c  T*  dn,  the  thermal  potential 

D  ^  dn,  the  dissipation  function 


(2-4) 


and 


T(t) 


,  the  thermal  force 


n=o 


which  must  satisfy  the  Lagrangian 


8 V  ^  ap 

8  qi  8qj 


Q 


(2-5) 


Solution  of  the  above  equation  then  requires 

•  specifying  the  heat  input  at  the  surface,  thus  determining 
the  thermal  force,  Q,  and 

•  specifying  the  spatial  temperature,  T,  so  that  the  function 
H  and  hence  the  functions  V  and  D  may  be  determined. 

The  solution ofEquation(2-5)willyieldthe  penetration  depth,  qj,  as  a  function 
of  time  and  hence  the  spatial  distribution  as  a  function  of  time. 

Because  of  the  Lagrangian  Equation  (2-5)^  the  form  of  the  temperature  distri¬ 
bution  specified  through  the  body  must  be  that  of  a  potential  function. 


2.  2.  2  External  Step  Function  Temperature  Input 


Applying  the  results  of  the  previous  section  to  a  circular  port  grain  with  an 
insulated  port  and  a  step  function  temperature  input  on  the  outside  requires 
specifying  the  spatial  distribution  of  T,  temperature. 

Referring  to  Figure  2-1.  the  geometrical  coordinates  are 

b  =  outer  radius  of  grain 

a  =  inner  radius 

r  =  radius  at  any  point. 
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Two  time  regimes  are  defined: 

•  The  time,  t,  range  o  <  t  <  t^#  where  t^  is  the  time  necessary 
for  the  heat  front  (defThed'as  the  distance  qi  from  the  outside 
to  reach  the  inner  port 

•  The  time  range  £  *  5  “  after  the  heat  front  has  reached 
the  inner  port 

The  input  function  for  the  step  temperature  is 


T(t) 


U(t) 


(2-6) 


where  U(t)  is  the  unit  step  function. 

The  spatial  temperature  is  assumed  in  the  quadratic  form  for  the  time 
range  until  the  heat  front  reaches  the  inner  port  as 


T{r) 


1  - 


O 


b  -  r  <  qi  and  o  <  t  <  t^ 
b  -  r  >  qi  and  t^  <  t  <  o 


(2-7) 


Rewriting  Equations  (2-1)  through  (2-5)  in  polar  coordinates,  defining 
Equation  (2-3)  for  H  in  terms  of  Equation  (2-6)  for  T^^  at  r  =  b,  defining 
V  and  D  in  terms  of  T(r)  and  H,  and  solving  Equation  (2-5)  yields  the 
relation 


q  q  -  2bq  q  =  -10Kb 
1  1  1  1 


(2-8) 


where 


K  =  ,  thermal  diffusivity,  ft^/sec 


for  the  penetration  depth,  q  ,  in  terms  of  k,  b,  and  t  (Figure  2-2), 
Solving  the  differential  Equation  (2-8)  for  q^  then  yields  the  cubic 


-  bq^  =  -10  K  bt 

r  1 


(2-9) 


whose  one  real  root  of  significance  is 


q^  =  b[2  cos  (iy-^)  +  1]  =  b  y  1  1.03  71^  =  Y 


(2-10) 
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where 

cos  6  =  1  -  15  T)* 

and 

“H*  =  ^  .  thermal  diffusion  parameter 
Rewriting  Equation  (2-10)  and  Equation  (2-12)  as 


(2-11) 


(2-12) 


10  / - 

-  "J-V  Kt 

it  can  be  seen  that  the  time,  ta.  for  the  heat  to  reach  the  inner  port,  q  =b-a 
is, letting  X  =  b/a  * 


,  0.09b*  /  X-  1,2 

or  somewhat  more  accurately 


=  ^( 


b*  /  2X*  -  3X*+  1 


) 


In  the  limit  as  a->-  o  and  X- 
t 


for  a  solid  cylinder 


(2 -13a) 


(2 -13b) 


and  for  a  thin  cylinder  as  a->-  b  and  X  -v  1 


O  . 


It  is  more  convenient  to  work  with  the  thermal  diffusion  parameters  instead 
of  the  time,  since  the  solution  is  then  nondimensionalized  to  facilitate  its 
extension  to  different  cylinder  sizes  and  materials.  From  Equation  (2 -13b) 
then 


^2  =  1  -  3X^  +  2X^ 

a  ioX’  ’ 

As  X  approaches  1.0,  1)1  approaches  zero, 
an  infinitely  thin  cylinder  instantaneously. 


(2-14) 

That  is,  the  heat  front  traverses 
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Once  the  temperature  field  reaches  the  inner  port  a  new  function  must  be 
selected  in  lieu  of  Equation  (2-7).The  new  function  must  meet  the  following 
requirements: 

•  At  r=a  the  temperature  must  gradually  increase  until  the 
temperature  of  the  entire  cylinder  approaches  a  constant 
value  (Figure  2-3). 

•  Since  the  inner  port  was  assumed  insulated,  the  derivative 
of  the  temperature  with  respect  to  r  must  be  zero  at  r=a. 

For  all  t  >  ta  or  >  T)j^  the  following  thermal  field  was  assumed: 


T  =  T, 


Cl 

•^o 


3  [1 


(b- 


+  T, 


(2-15) 


where  Ta  is  defined  in  Figure  2-3. 


Setting  r=a,  it  is  clear  that  Ta  represents  the  temperature  at  the  inner  port, 
and,  hence,  is  a  function  of  time.  The  procedure  for  finding  this  functional 
relationship  is  similar  to  that  used  for  qj,  yielding 


Ta  =  To  [1  -  ] 

where 

10X^[3X+  1]  1)1  C3X.f  1]  C2\+  1] 

■  (X-l)*  [5X+1]  "  3X  [5X+11] 


and  where 


(2-16) 


(2-17) 


An  inspection  of  Equation  (2-16)  reveals  that  when 


7)2  =  ;  Ti"  =  1 

then 


Ta  =  0 
when 
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Figure  2-3  Change  of  Inner  Port  Temperature  with  Time 
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Since,  as  stated  above,  represents  the  temperature  of  the  inner  port, 
the  relation  for  is  physically  reasonable.  The  results  for  a  cylinder 
with  b/a  equal  to  2  are  shown  in  Figure  2-4.  Equations  (2-7)and  (2-15)  com¬ 
pletely  describe  the  temperature  distributions  in  a  hollow  cylinder  with  an 
insulated  inner  port  and  a  case -temperature  variation  consisting  of  a  step 
change  T  .  The  values  of  t^,  qj  T^,  and  7  can  be  determined  from  Equations 
(2-5),  (2-13),  (2-16),  and  (2-17)  for  g  i  ven  values  of  X  and  k  .  The  para¬ 
meter  X  depends,  of  course,  on  the  geometry  and  the  thermal  diffusivity  k 
on  the  type  of  propellant.  Plots  of  1^,  7 ,  k  ,  and  T^  are  the  only  required 
input  parameters  necessary  for  the  thermal  field  determination  in  any 
hollow  cylinder  with  the  aforementioned  boundary  conditions. 


2.  2.  3  Arbitrary  External  Temperature  Variation 


In  determining  thermal  fields  for  more  general  types  of  case -temperature 
variations,  the  following  approach  was  used: 

•  Representation  of  an  arbitrary  surface  temperature  variation 
by  a  series  of  step  changes.  This  allows  the  direct  use  of  the 
step  solution  described  in  Section  2.  2.  2.  The  solution  for  a 
step  surface  temperature  is  known;  hence,  if  a  varying  surface 
temperature  is  degenerated  into  discrete  steps,  the  resulting 
thermal  field  can  be  accurately  approximated' by  correctly 
adding  the  appropriate  step  solutions. 

The  use  of  this  method  is  illustrated  in  the  following  example.  Consider 
as  before  a  hollow  cylinder,  initially  at  thermal  equilibrium,  subject  to  a 
step  exterior  surface  temperature  change  at  time  zero.  At  the  time  when 
the  resulting  heat  front  arrives  at  the  inside  of  the  cylinder,  a  step  exterior 
temperature  change  opposite  in  sign  to  the  first  condition  is  imposed.  The 
temperature  input  and  resulting  thermal  fields  in  the  cylinder,  obtained  by 
superposition  of  the  two-step  input,  are  shown  in  Figure  2-5.  Figure  2-5 
indicates  the  complexity  of  interaction  of  the  two  heat  inputs  leading  to  a 
transient  condition  where  heating  and  cooling  proceed  simultaneously  in 
different  parts  of  the  cylinder.  It  is  this  interaction  which  makes  cyclic 
inputs  an  interesting  problem  for  viscoelastic  materials. The  analysis  for 
this  problem  is  given  in  Appendix  A. 

The  method  illustrated  by  the  hollow  cylinder  example  discussed  above  can 
be  systematically  extended  to  a  sequence  of  step  inputs  without  introducing 
any  restrictive  complexity,  and  was  used  for  the  initial  evaluation  of  the 
viscoelastic  thermal  stress  accumulation  problem. 


2-  1 1 


UOCKHCEO  PROPULSION  COMPANY 


2-12 


LOCKHEED  PROPULSION  COMPANV 


Figure  2-4  Time -Temperature  Distribution  in  a  Hollow  Cylinder  (port  insulated)  with 
Zero  Initial  Temperature  and  a  Suddenly  Applied  Surface  Temperature,  T, 
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2.  3  PROPELLANT  THERMAL  PROPERTIES 
2.  3.  1  Conductivity,  Specific  Heat  and  Diffusivity 


Before  numbers  can  be  extracted  from  the  thermal  field  solutions,  the  thermal 
properties  must  be  obtained.  The  pertinent  properties  are: 

c  =  thermal  conductivity 

K  =  thermal  diffusivity 

k  =  specific  heat 


where 


K  =  k/ pc  ,  p  =  density 

K,  k,  and  c  are  dependent  on  temperature.  However,  the  inclusion  of  such 
a  dependence  in  the  heat  conduction  equation  would  result  in  a  nonlinear 
equation  whose  solution  would  be  too  complicated  to  be  of  value  for  the  pre¬ 
sent  study.  The  experimental  effort  at  LPC  was  directed  at  obtaining  engin¬ 
eering  values  rather  than  scientific  values  for  the  thermal  properties.  It 
was  felt  that  this  approach  would  expedite  the  extraction  of  sufheienily 
accurate  values  for  the  properties  for  use  in  the  thermal  field  calculations. 
The  experimental  values  obtained  for  these  constants,  when  used  in  the 
theoretical  thermal  field  solutions,  did  in  fact  result  in  thermal  field  pre¬ 
dictions  which  agreed  accurately  with  the  thermal  field  tests  in  the  tempera¬ 
ture  range  of  interest  (70®  to  -100®F). 


The  experimental  values  for  c,  the  thermal  conductivity,  were  obtained 
using  a  standard  Cenco-Fitch  therm"!  conductivity  apparatus  (Figure  2-6). 
Tests  of  Polycarbutene  R  propellant  gave  an  average  value  for  c  of 


c 


9.34  X  lO"* 


Cal -Cm 
Sec-C|n-®c 


Measurements  of  the  specific  heat  were  performed  using  a  standard  labora¬ 
tory  calorimeter.  The  specimens  were  coated  with  silicone -base  paint  to 
prevent  absorption  of  the  conducting  medium  by  the  propellant  sample.  An 
average  value  for  the  specific  heat  as  obtained  from  replicate  tests  was 

k  =  0.26  BTU/lb/®F 


The  resulting  thermal  diffusivity  was  calculated  to  be: 
K  =  134.1  X  10‘‘ ftVmin. 
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Figure  2-6  Cenco-Fitch  Thermal  Conductivity  Apparatus 
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Z.  3.  2  Linear  Expansion 


Experimental  apparatus  for  measuring  the  coefficient  of  thermal  expansion 
consisted  of  a  simple  quartz -tube  linear  dilatometer  employing  electronic 
means  for  length  measurement  (Figure  2-7). 

The  linear  expansion  thermal  coefficient  of  Polycarbutene  propellant  in  the 
7C®  to  -100®F  range  was  5.6  ±  0.2  in./4n./®F. 


2.  3.  3  Heat  Transfer  Coefficients 


The  determination  of  thermal  fields  for  a  motor  placed  suddenly  into  a  con¬ 
ditioning  box  involves  heat  transfer  through  an  air  boundary  layer  and  heat 
conduction  through  a  composite  cylinder  composed  of  steel,  rubber  insula¬ 
tion,  and  a  hollow  grain.  The  heat  transfer  coefficient  is  a  function  of 
geometry,  air  temperature  and  air  circulation.  Each  motor  would  have  a 
different  heat  transfer  coefficient;  however,  for  practical  considerations, 
a  reasonable  average  value  for  the  heat  transfer  coefficient  could  be  used 
for  different  motors. 

Two  tests  were  performed  at  LPC  to  determine  a  suitable  value  for  the  heat 
transfer  coefficient  for  a  4  by  10-inch,  ‘/t-inch  thick,  steel  case.  The  inter¬ 
ior  of  the  case  was  stuffed  with  insulating  material  and  a  thermocouple  was 
embedded  in  the  case.  The  resulting  case  temperature  versus  time  was 
plotted.  An  approximate  value  for  the  heat  transfer  coefficient  was  then 
calculated  based  on  the  following  assumptions: 

•  The  temperature  gradient  through  the  case  is  small  so  that 
the  case  has  a  uniform  but  time -dependent  temperature  and 
no  heat  is  lost  through  conduction  into  the  interior. 

•  The  heat  input  to  the  case  is  directly  proportional  to  the 
difference  in  the  ambient  and  case  temperature. 

These  assumptions  lead  to  the  following  relations: 

q  =  h{To  -  Tc)  (2-18) 

where 

q  =  heat  input  in  BTU/hr 
h  =  heat  transfer  coefficient  in  BTU/hr  -  ®F 
Tq  =  temperature  of  ambient  air  ®F 
Tc  =  temperature  of  case  ®F 
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and 

^  (2-19) 

where 

M  =  total  mass 
Cp  =  specific  heat  of  case 

Substituting  Equation  (2-18)  into  Equation  (2-19)  results  in 

h(To  -  Tc)  =  MCp-^^'^o  ■ 
hence 

ijr  _  I’p-I'c  _  ^-ht/MCp 
^o 

The  relations  indicate  that  a  plot  of  T"  should  be  a  straight  line  on  semi -log 
paper  and  the  value  of  h  can  be  determined  from  the  slope  of  this  line.  The 
plot  of  the  experimental  data  obtained  from  tests  on  the  steel  case  mentioned 
above  did  in  fact  prove  to  be  a  straight  line  on  semi -log  paper  (Figure  2-8). 
The  value  for  h  as  obtained  from  these  tests  was 

h  =  5.55®TU/“F-hr 


(2-20) 

(2-21) 
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Figure  2-8  Heat  Transfer  Coefficient  Test 
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2.4  EXPERIMENTAL  TEMPERATURE  DISTRIBUTIONS  IN 
CIRCULAR  PROPELLANT  GRAINS 


Experimental  tests  were  performed  at  LPC  to  verify<the  approximate  varia¬ 
tional  step  solutions  and  other  theoretical  solutions  based  on  them,  as  well 
as  the  accuracy  of  the  thermal  properties. 

Thermocouples  were  cast  in  situ  in  6-inch  diameter  cylindrical  grains  of 
Polycarbutene  R  propellant  (Figures  2-9  and  2-l0).  The  grains  were  solid 
cylinders  with  heavily  insulated  ends.  During  the  tests  the  thermocouples 
were  monitored  to  obtain  the  time-temperature  variation  at  various  radii  in 
the  grains. 


2.4.1  External  Step  Temperature  Variation 


The  step  surface  temperature  variation  test  was  performed  by  immersion  of 
a  laquer  coated  grain,  initially  at  thermal  equilibrium  at  70°F,  in  a  liquid 
bath  at  20®F.  The  resulting  temperature  versus  dimensionless  time  as 
measured  at  three  radial  positions  in  the  grain  are  compared  with  the  exact 
and  approximate  solutions  in  Figures  2- 1 1 ,  2-12,  and  2-13.  The  experimen¬ 
tal  values  agree  within  5  or  6  percent  with  those  predicted  by  the  approximate 
technique. 


2.4.2  Arbitrary  External  Temperature  Variation 


Experiments  in  which  sequential  step  and  continuously  varying  thermal  inputs 
were  imposed  on  a  solid  cylinder  were  performed  for  evaluation  of  the  super¬ 
position  me  th  o  d  o  f  s  p  e  ci  f  yin  g  th  e  r  mal  fields  f  o  r  time  varying 
thermal  inputs.  In  the  sequential  or  double  step  experiment,  a  coated  solid 
cylinder  was  subjected  to  a  step  surface  temperature  variation  from  equilib- 
rivim  at  70°F  to  20^F  by  immersion  in  fluid  at  the  latter  temperature.  After 
conditioning  at  20°F  for  45  minutes,  the  cylinder  was  immersed  in  a  second 
liquid  bath  at  120'’f.  The  temperature  variation  with  time  was  measured  at 
different  radial  positions  in  the  cylinder.  In  Figures  2-14,  2-15  and  2-16, 
data  for  the  test  are  compared  to  predictions  based  on  the  step  solution  super¬ 
position  approach  previously  described.  The  agreement  between  the  predicted 
and  measured  temperature  time  variation  was  within  5  percent. 

Evaluation  was  made  of  the  effects  of  a  continuously  varying  surface  tempera¬ 
ture  on  a  solid  grain  by  a  test  in  which  a  grain  was  first  subjected  to  a  step 
decrease  in  surface  temperature  to  SS'^F  by  liquid  immersion.  Subsequently, 
the  surface  temperature  of  the  grain  was  warmed  to  120*’F  in  the  liquid  bath 
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Figure  2-9  Sealrite  Casting  Assembly  (side  view) 
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Figure  2-10  Sealrite  Casting  Assembly  (top  view) 
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Figure  2-11  Temperature  versus  Time  in  a  Solid  Cylinder  Subjected 
to  a  Sudden  Temperature  Change  at  r  =  0.  45  b 
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Figure  2-12  Temperature  versus  Time  in  a  Solid  Cylinder  Subjected 
to  a  Sudden  Temperature  Change  at  r  =  0.45  b 


Figure  2-13  Temperature  versus  Time  in  a  Solid  Cylinder  Subjected  to 
a  Sudden  Temperature  Change  at  r  =  0.  896 
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Figure  2-15  Temperature  versus  Time,  r  =  0.466,  Sequential 
Step  Surface  Temperature  Variation 
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igure  2-16  Temperature  versus  Time  at  r  =  0.896,  Sequential 
Step  Surface  Temperature  Variation 
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over  a  period  of  155  minutes  and  held  constant  at  120**F  (Figure  2-17).  The 
temperature -time  variation  at  two  radii  in  the  grain  is  compared  with  that 
predicted  by  the  superposition  integral  approach  previously  delineated  in 
Figures  2-18  and  2-19.  The  observed  agreement  was  again  within  5  percent. 
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2.5  DISCUSSION 


While  the  experimental  thermal  field  data  for  propellant  grains  examined  in 
this  program  are  not  sufficient  for  a  realistic  error  analysis,  the  results 
show  that  good  engineering  accuracy  can  be  obtained  by  using  the  approximate 
energy  techniques.  The  problems  considered,  theoretically  and  experimen¬ 
tally,  treated  only  the  condition  where  a  uniform  surface  temperature  of  an 
ixisymmetric  body  was  specified.  The  prediction  of  thermal  fields  in  a 
rocket  grain  heated  or  cooled  with  heat  transfer,  through  a  gaseous  external 
boundary  layer  is  subject  to  considerably  larger  error  because  of  the  un¬ 
certainty  of  the  heat  transfer  coefficients.  A  second  source  of  error  in  cal¬ 
culation  of  thermal  fields  in  rocket  grains  in  heat  transfer  through  the  motor 
case,  liner  and  insulation.  It  is  expedient  and,  indeed,  sufficiently  accurate 
for  engineering  purposes  to  ignore  thin  liners  and  insulation  and  to  assume 
no  thermal  gradient  through  the  case.  A  comparison  of  predicted  and 
observed  thermal  fields  for  this  situation  are  shown  in  Section  5.2. 

The  problem  of  axi symmetric  thermal  fields  in  star  perforated  propellant 
grains  is  accessible  to  solution  by  using  generalized  thermal  field  analysis 
techniques  based  on  digital  computer  methods.  The  associated  transient 
viscoelastic  stresses  for  the  star  grain  cannot,  however,  be  rigorously  spec¬ 
ified  except  by  approximation  of  the  grain  as  a  circularly  perforated  one  with 
stress  concentrations.  The  generalized  problem  of  arbitrary  non-socisymmet- 
ric  thermal  inputs  to  a  rocket  grain  is  as  yet  an  unresolved  problem  -  -as  is 
the  associated  viscoelastic  stress  analysis  problem. 

The  present  problem  of  axisymmetric  thermal  fields  in  an  infinite  length 
circ^arly  perforated  grain  has  been  resolved  to  an  acceptable  degree  to 
engineering  practice  in  a  theoretical  framework  uncomplicated  enough  to 
permit  introduction  of  the  thermal  field  data  into  the  transient  viscoelastic 
stress  analysis  framework. 
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Section  3 

PROPELLANT  VISCOELASTIC  CHARACTERIZATION 


3.  1  BACKGROUND 


Solution  of  the  transient  thermal  stress  problem  requires  the  assumption  of 
thermorheologically  simple  behavior  for  the  propellant,  which  is  treated  as 
a  linear  viscoelastic  material.  Thermorheological  simplicity  stipulates  a 
relationship  between  the  effects  of  time  and  temperature  whereby  a  change  in 
temperature  produces  a  change  in  a  given  stress-strain  relationship  identical 
to  that  obtained  if  the  time  scale  of  observation  were  to  be  expanded  or  con¬ 
tracted  without  alteration  of  the  temperature.  In  view  of  the  importance  of 
the  thermorheological  viscoelastic  relationships  in  the  transient  thermal 
stress  problem,  the  temperature -time  reduced  variable  approach  was  exam¬ 
ined  in  considerable  detail.  The  development  of  reduced  variable  mathemat¬ 
ical  treatments  for  operational  and  hereditary  forms  of  the  stress-strain 
law  in  uniaxial  tension  is  presented  in  detail  in  Section  3.  2.  Experimental 
data  for  the  Polycarbutene  R  Propellant  studied  in  this  program  are  included 
in  Section  3.3. 
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3.  2  REDUCED  VARIABLE  CONCEPTS 


The  concept  of  reduced  variables  or  viscoelastic  corresponding  states  pro¬ 
vides  a  valuable  simplification  in  the  mathematical  description  of  the  mechan¬ 
ical  behavior  of  linear  viscoelastic  materials  by  separating  the  principle 
variables  time  and  temperature.  In  reduced  variable  representation  the 
properties  of  a  polymer  are  represented  in  terms  of  single  functions  of  time 
and  temperature  whose  forms  can  be  determined  experimentally  even  when 
the  function  representing  temperature  effects  is  not  conveniently  represented 
in  analytic  form  (Ferry,  1961). 

In  its  least  restrictive  sense,  the  concept  of  viscoelastic  corresponding  states 
can  be  defined  as  follows.  Mechanical  behavior  of  a  linear  viscoelastic  mate¬ 
rial  at  different  equilibrium  temperatures  can  be  related  through  treatment  of 
behavior  at  the  various  temperatures  in  terms  of  a  reduction  parameter  which 
stipulates  a  numerical  relationship  between  the  relaxation  and  revardation 
spectra  for  the  material  as  a  function  of  temperature.  The  reduction  para¬ 
meter  specifies,  for  example,  the  ratio  of  relaxation  times  at  one  tempera¬ 
ture  to  the  corresponding  relaxation  times  at  another  temperature. 

Through  time -temperature  superposition,  hence,  data  for  various  tempera¬ 
tures  can  be  translated  in  stress,  strain,  time  space  to  superimpose  and 
form  a  single  or  reduced  representation  of  stress-time  or  strain-time  behav¬ 
ior  at  an  arbitrary  reference  temperature.  Reduced  or  master  physical 
property  plots  will  in  general  display  the  viscoelastic  behavior  of  the  material 
over  a  wider  span  of  time  than  can  be  obtained  in  experiments  possible  at  a 
single  temperature.  In  tests  of  a  viscoelastic  material  at  one  temperature, 
observation  of  effects  associated  with  a  very  short  relaxation  time  is  ham¬ 
pered  by  the  experimental  necessity  for  causing  finite  stress  or  strain  vari¬ 
ations  in  a  time  interval  of  the  same  order  as  that  of  the  relaxation  or 
retardation  time  itself.  The  accurate  elucidation  of  relaxation  times  of  below 
10"^  seconds  at  a  given  temperature  is  an  experimental  problem  of  consider¬ 
able  difficulty.  However,  a  reduction  in  the  test  temperature  of  a  practical 
amount  will  in  effect  slow  the  short  relaxation  times  so  as  to  bring  them  into 
the  easily  observable  experimental  range.  Hence,  the  reduction  in  tempera¬ 
ture  causes  a  "magnification"  of  the  effects  of  the  short  relaxation  times  both 
in  terms  of  their  time  duration  and  their  contribution  to  stress  or  strain 
variation  on  the  reduced  time  scale.  The  converse  effects  brought  about  by 
increasing  temperature  are  observed  also;  relaxation  or  retardation  times 
too  slow  to  be  observed  at  one  temperature  can  be  speeded  by  increases  in 
temperature. 
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For  materials  that  are  thermorheologically  simple  (i.  e.  ,  those  that  can  be 
treated  accurately  in  reduced  variable  form)  the  concept  of  viscoelastic  cor¬ 
responding  states  is  an  empirical  adjunct  of  vital  importance.  Experimental 
data  thai  c^n  be  resolved  accurately  in  a  moderate  time  span  only  can  be 
obtained  at  various  temperatures  and  reduced  to  provide  data  for  a  broad 
spectrum  determination  of  the  viscoelastic  properties  of  a  material.  Broad 
spectrum  viscoelastic  material  property  data  are  essential  for  stress  analy¬ 
sis  of  rocket  grains,  both  for  consideration  of  dynamic  loading  and  for  the 
transient  thermal  stress  problem. 

The  following  presents  an  approach  to  the  mathematical  .treatment  of  the 
reduced  variable  concept.  The  mathematical  treatment  is  simplified  for 
consideration  of  the  reduced  variable  method  in  power  law  and  linear  differen¬ 
tial  operator  viscoelastic  material  specifications,  and  for  development  of 
methods  to  accurately  determine  the  temperature -time  reduction  variable 
through  experimentation. 


3.  2.  1  Basic  Reduced  Variable  Relationships 


The  reduced  variable  treatment  of  time -temperature  viscoelastic  properties 
takes  into  account  the  variation  of  the  viscous  and  elastic  aspects  of  polymer 
behavior  with  absolute  temperature.  The  terms  viscous  and  elastic  are  used 
herein  in  a  descriptive  sense  rather  than  in  an  explicitly  interpretative  sense 
with  respect  to  polymer  molecular  mechanics.  The  temperature -time  viscous 
behavior  is  related  through  the  shift  parameter  a.j.  by  noting  that  the  recipro¬ 
cal  relaxation  times  Q'^(T)  for  a  material  at  absolute  temperature  T  are 
related  to  the  relaxation  times  ®i{Tg)  at  an  arbitrary  reference  temperature 
Tg  by 

Q'i(T)  =  aT/a'.(Ts)  (3-1) 

where  a-p  by  definition  equals  unity  at  the  reference  temperature.  The  elas¬ 
tic  component  of  behavior  is  related  to  temperature  through  noting  that  the 
elastic  component  Ej(T)  at  temperature  T  is  related  to  the  elastic  component 
Ej(Tg)  at  the  reference  temperature  Tg  by 

Ei(T)  =  ^  Ei(Tg)  (3-2) 

^  3 


This  latter  observation  (Equation  3-2)  arises  theoretically  from  considera¬ 
tion  of  an  equilibrium  molecular  rubber  network  as  a  network  of  entropy 
springs  (Treloar,  1958).  The  validity  of  the  relationship  expressed  by 
Equation  3-2  has  been  demonstrated  experimentally  in  tests  on  fully  relaxed 
rubbers.  The  validity  of  the  shift  parameter  a-p  has  been  demonstrated 
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empirically.  Theoretical  rationales  for  interpretation  of  the  molecular 
origin  and  significance  of  aj  have  been  advanced,  but  their  consideration  is 
not  germane  to  this  discussion.  Theoretical  rationale  has  also  predicted  the 
existence  of  second  order  time -temperature  reduction  parameters,  such  as 
a  density  ratio  for  various  temperatures.  The  second  order  reduction  para¬ 
meters  will  be  neglected  in  the  subsequent  discussion  since  they  are  clearly 
below  the  level  of  numerical  significance  in  engineering  treatments  of  prob¬ 
lems  pertinent  to  this  contract. 

Separation  of  that  portion  of  the  elastic  response  associated  with  glassy 
behavior  as  a  temperature  independent  phenomenon,  as  has  been  suggested 
(Ferry,  1961),  will  be  ignored  also.  While  there  are  theoretical  approaches 
that  suggest  this  is  advisable,  in  practice  the  separation  of  the  effect  is  com¬ 
plicated  and  the  associated  effects  are  trivial  in  engineering  calculations. 


3.  2.  2  Reduced  Variable  Treatment  for  Linear  Viscoelastic  Materials 


A  viscoelastic  material,  linear  or  not,  cannot  be  considered  to  be  a  structure 
of  springs  and  dashpots  in  the  literal  sense.  However,  the  mathematics 
describing  linear  viscoelastic  physical  behavior  also  describe  the  mechanical 
behavior  of  a  network  of  Hookean  springs  interconnected  with  Newtonian 
dashpots.  Hence,  because  of  the  inherent  analogy  between  linear  viscoelastic 
behavior  in  the  real  sense  and  the  behavior  of  a  linear  network  composed  of 
elastic  springs  and  viscous  dashpots,  it  is  both  expedient  and,  with  some 
restrictions,  mathematically  valid  to  rationalize  linear  viscoelastic  behavior 
in  terms  of  a  spring  and  dashpot  assemblage. 

Consider,  for  example,  the  model  analog  illustrated  in  Figure  3-1, 


Figure  3-1  Spring  and  Dashpot  Analog 
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where  Ej,  E2  are  the  stress-strain  constants  of  Hookean  springs  and  1)2  is  a 
Newtonian  viscosity.  For  convenience,  force -deflection  relationships  will  be 
considered  in  terms  of  the  normalized  quantities  stress  and  strain.  In  gen¬ 
eral  operational  form,  the  stress  strain  law  for  the  analog  is 

where  CT(t)  and  e  (t)  are  the  time  dependent  stress  and  strain,  respectively. 
The  terms 


111 

El  +  E2 


1_ 

O' 


and 


Ez 


(3-4) 


(3-5) 


are  the  relaxation  and  retardation  times 
empirical  reduced  variable  approach,  it 
as 

d  Tr(Ei+E2) 

"51^  TI2  At 


a(T)  =  Tr 


for  the  analog.  Consistent  with  the 
is  appropriate  to  write  Equation  (3-3) 


E, 


Tr  Ez 

T)  2  At 


e(T)  (3-6) 


where  t  is  reduced  time,  Tj.  is  the  ratio  of  absolute  temperature  T  to  the 
reference  temperature  Tg  and  At  is  defined  as 


(3-7) 


The  differentiation  between  At  as  used  here,  and  aTi  as  developed  historic¬ 
ally,  is  made  primarily  to  indicate  a  separation  between  absolute  temperature 
effects  on  the  "springs"  and  "viscosity"  changes  with  temperature  in  the 
"dashpots". Solution  of  Equation  (3-6)  for  a  step  function  strain  input  condition 
at  an  equilibrium  temperature  T  yields,  for  reduced  time: 


a(T)  =  Tj.  E, 


e  -  (e  -ff)e’®* 
a 


Jt 

AT 


(3-8) 
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In  proceeding  to  consideration  of  the  generalized  solution  of  more  pertinent 
engineering  interest,  it  is  convenient  to  consider  the  model  made  up  in  Fig¬ 
ure  3-2  below,  of  parallel  assemblages  of  the  model  shown  in  Figure  3-1. 


Figure  3-2  Parallel  Assemblages  of  Model  Analog 


The  generalized  operator  equation  for  this  model  is 


3"  a^-‘ 

—  +^n-i  TTTTT  +---+ao 


at 


St* 


a(t) 


n-i 


+  b 


at' 


n-i  -5rr  «(t) 


at 


1 


(3-9) 


where  the  coefficients  of  the  operator  a^,  bj  describe  the  specific  behavior  of 
the  generalized  model.  In  terms  of  the  model,  the  coefficient  bj^  is  identified 
as 


n 

=  2 
i=i 


'1.  i 


(3-10) 
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Equation  (3-9)  can  be  solved  for  either  analytic  stress  or  strain  input 
functions.  The  input  functions  of  general  interest  are: 

Stress  Relaxation:  e  (t)  =  Cq  U  (t) 

where  U  (t)  is  the  Heaviside  unit  step  at  t  =  o  and  Cq  the 
constant  strain  level; 


Creep:  Oq  U  (t) 

where  Oq  is  the  constant  stress  level; 
Constant  Strain  Rate:  c  (t)  =  Rt  U  (t) 


where  R  is  the  strain  rate; 

Ramp  Strain  Relaxation;  c  (t)  =  R  t  U{t)  -  R  (t  -  ti)  U  (t  -  ti) 

where  U  (t  -  ti)  is  the  Heaviside  unit  step  at  t  =  tj.  Physically,  the  ramp 
stress  relaxation  input  is  a  constant  strain  rate  ramp  to  t  =  ti  and  a 
constant  strain  level  e(ti)  =  Rti  for  t  >  ti.  The  solutions  for  these  input 
conditions  are  given  below  for  the  reduced  variable  interpretation: 


Stress  Relaxation 


^  ,  1 


Creep 

e(T)  =  ;^ 


B 


+.T  B,  e-^i 

n  1=1  1 

'i>n  S; 

1=1 

(S  +  ak) 


b.n  S  n  ,  (S  +  0k) 

k=i 


(S  +  0i) 


s  =  -ei 

3i/^o 


(3-13) 


(3-14) 


(3-15) 
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Constant  Strain  Rate 

"(^)=RTy  S  E,.i 
i=i 


Ramp  Stress  Relaxation 


Interpreting  these  relationships  physically,  it  is  seen  that 
stress  relaxation  data  for  various  temperatures  plotted  as 


(3-16) 


q(t) 

fioTr 

T 


vs ,  t 


(3-18) 


and  creep  data  plotted  as 


r.(t) 


vs .  t 


(3-19) 


Can  be  shifted  along  the  time  scale  to  determine  the  shift 
factor  a^ 


At 

Tr 


(3-20) 
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with  reference  to  an  arbitrary  reference  temperature  Tg  at 
which  a^p  is  defined  as  unity.  Constant  strain  rate  data  plotted 
as 


ojt^ 

RTj. 


Aq  t 


vs.  t , 


T 


(3-21) 


where  Ag  is  the  rubbery  equilibrium  modulus,  can  be  shifted 
at  45°  ( ap  vertically  and  ap  horizontally)  to  obtain  the  shift 
variable. 


In  reduction  of  the  ramp  stress  relaxation  data,  a  three  com¬ 
ponent  shift  is  necessary  to  superpose  data  for  t  >ti  because 
of  the  term 


(  1  -  e 


Oi 


(3-22) 


in  Equation  (3-17).  Rigorously,  the  shift  factor  cannot  be 
found  from  graphical  plots  of  the  ramp  stress  relaxation  data. 
On  a  practical  basis,  however,  the  solution  for  Equation (3- 17) 
degenerates  to  that  for  Equation  (3-13)  for  t  >  11  tj,  to  within 
a  fraction  of  a  percent  accuracy,  the  approximation  becoming 
more  accurate  as  t  increases  beyond  11  ti.  Thus,  accurate 
values  for  the  shift  factor  from  ramp  stress  relaxation  data 
may  be  obtained  by  simple  graphical  interpretation  of  the  data 
for  t  >  11  t,  in  terms  of  the  step  function  stress  relaxation 
solution. 

The  reduced  variable  plots  for  stress  relaxation  and  creep. 


^  (^)  vs.  Va 

^oTr 


and 


vs. 


(3-23) 


(3-24) 


yield  the  master  relaxation  modulus  and  creep  compliance 
curves  for  the  material  at  the  reference  temperature  Tg . 
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The  reduced  variable  plot  for  constant  strain  rate  data 


ax  vs.  ‘/a^  (3-25) 

yields  the  master  stress- strain  curve  for  unit  strain  rate  for  the 
material  at  the  reference  temperature  Tg  .  The  time  derivative 
of  this  curve  is  equal  to  the  relaxation  modulus  at  Tg  .  In  each 
of  the  above  cases,  the  master  curves  for  the  moduli  or  com¬ 
pliance  together  with  the  time-temperature  shift  factor  awfully 
specify  the  material  properties  for  the  temperature  and  reduced 
time  domain  considered. 

Master  curve  data  at  the  reference  temperature  are  converted 
to  master  curve  data  at  another  temperature  T^^  by  "unreducing” 
the  data  in  accord  with  the  ratio  T^/Xg  the  value  of  a  ^  for 

the  temperature  T^^  . 


3.2.3  Power  Law  Approximations  to  the  Viscoelastic  Stress-Strain  Law 


The  heriditary  integral  form  of  the  stress- strain  law 


a(t)  = 


I 


H  {  %  )  a 


+  a) 


d  (  ya  )  _d_ 

"dt 


p  (t) 


(3-26) 


where 

Eg  =  equilibrium  modulus 
H(yQ')  =  relaxation  distribution 


is  an  alternate  form  equally  acceptable  mathematically  to  the  operational 
stress-strain  law  (Equation  3-9).  Indeed,  the  limit  consideration  of  the 
hereditary  form  by  a  power  law  approximation  is  of  particular  interest 
and  was  used  as  the  material  properties  statement  in  the  solution  of  the 
transient  thermal  stress  analysis  problems. 
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For  purposes  of  discussion,  letting  the  term  Eg  of  Equation  (3-Z6)  be  zero 
and  the  expression  under  the  integral  be  infinity  for  t  =  o,  then  the  power  law 
form  for  the  relaxation  modulus  Er  (t)  becomes 


E  r  (t)  =  r  Ki 


1  ”i 

(t)  • 


(3-27) 


For  the  class  of  propellants  studied  in  this  investigation  it  was  found  that  the 
above  summation  was  given  accurately  by 


Ej.  (t)  =  K  t  ^  (3-28) 

for  the  time  and  temperature  range  of  interest  in  the  transient  thermal  stress 
problem.  As  implied  by  Equation  (3-28),  the  stress  relaxation  modulus  was 
a  straight  line  on  a  plot  of  log  Ej.  (t)  versus  log  time. 

In  reduced  variable  form,  the  power  law  approximation  (Equation  3-28)  is, 
for  the  various  test  conditions  of  interest 


Stress  relaxation 


€ 


0 


n-1 


Creep 


<  (t) 

CTo 


(ik)  (^J 


(3-29) 


(3-30) 


Constant  strain  rate 


<j(t) 

R 


(3-31) 
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Inter  conversion  of  the  expressions  (3-29).  (3-30),  and  (3-31)  with  the  expres 
sions  (3-13),  (3-14),  and  (3-15)  over  the  range  where  the  power  law  approxi¬ 
mation  is  valid  can  be  done  for  known  or  assumed  values  of  or  8^  in 
Equations  (3-29),  (3-30)  and  (3-31).  As  a  more  direct  alternative  perhaps, 
a  single  set  of  data  can  be  fit  by  either  the  operational  solutions  or  by  the 
power  law  form  and.  in  the  range  of  applicability,  are  by  definition  equival¬ 
ent  statements.  Choice  of  the  method  for  properties  specification  is  a 
matter  of  convenience  in  application  of  the  mathematical  properties  state¬ 
ment  to  solution  of  the  problems  of  interest. 
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3.3  PROPELLANT  VISCOELASTIC  DATA 


The  uniaxial  viscoelastic  stress-strain  law  for  Polycarbutene  R  propellant 
was  obtained  from  interpretation  of  constant  strain  rate  data  for  the  propellant. 
Tab-end  tensile  (specimens  of  the  propellant  were  tested  in  replicate  at  tem¬ 
peratures  between  -100®  and  +  142®F.  The  data  were  superposed  for  a  ref¬ 
erence  temperature  of  70®F  according  to  Equation  (3-21)  to  obtain  the  master 
stress-strain  curve  for  a  strain  rate  of  0.005  minute  ■*  (Figure  3-3)  and  the 
associated  shift  factor  variation  with  temperature  (Figure  3-4).  The  strain 
levels  in  all  tests  were  confined  to  values  below  10  percent. 

The  power  law  approximation  for  the  stress  relaxation  modulus,  in  reduced 
time,  was 


Er  (t)  =  Tj.330  /  t_ 

rT/ 


(3-32) 


10®  >  t/ 

/a-p 


>  10 


and  the  operational  form,  as  expressed  by  the  form. 


II 

u 

n 

Tr  r  A. 
i=o 

e  -  O'  i  t/a,p 

i 

Ai 

Oj 

0 

63 

0 

1 

40 

10-3 

2 

66 

10-^ 

3 
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10-* 

4 
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1 

5 

290 

10 

6 

476 

10* 

7 
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10* 

8 

1,280 

10® 

9 

2,100 

10* 

10 

3,500 

10® 

11 

5,700 

lo’ 

12 

9,300 

10® 

13 

15,300 

lo’ 

14 

25,000 

10*® 

15 

41,000 

10** 

16 

67,000 

10** 
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In  consideration  of  the  nature  of  the  transient  thermal  stress  problem,  the 
introduction  of  the  absolute  temperature  ratio  in  the  properties  statement 
was  considered  to  be  of  trivial  importance  in  engineering  calculations.  Hence, 
the  material  properties  statement  as  entered  in  the  thermal  stress  calcula¬ 
tions  was 


Er  (t) 


10^  >  t  >  10’“ 


(3-34) 


and  the  associated  shift  factor  values  a'rp  measured  for  the  data  considering 
the  value  to  be  unity  for  all  temperatures.  Comparative  values  for  the 
shift  factors  as  derived  by  the  two  techniques  are  listed  in  Table  3-1. 
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Table  3-1 

TEMPERATURE  -  TIME  REDUCTION  PARAMETERS 


t("f)  ^T  ^T 


115 

2.95 

X 

10"^ 

4.31  X  10 

70 

1 

1 

20 

8.93 

X 

10* 

7.04  X  10 

-16 

1.04 

X 

10^ 

4.51  X  10 

-20 

2.70 

X 

10* 

1.14  X  10 

-50 

3.33 

X 

10* 

1.01  X  10 

-70 

2.03 

X 

10* 

4.85  X  10 

-80 

3.54 

X 

10* 

7.59  X  10 

100 

3.35 

X 

10** 

5.54  X  10 

3-17 


LOCKHKBD  PROPULSION  COMPANY 


578-F 


3.5  DISCUSSION 


The  preceding  discussions  have  not  lead  to  the  conclusion  that  reducedvari- 
able  treatment  of  propellant  properties  for  stress  analysis  problems  is  a 
completely  resolved  problem.  However,  in  view  of  the  successful  reduced 
variable  treatment  of  the  Polycarbutene  R  propellant  studied,  it  seems  clear 
that  within  engineering  accuracy  the  reduced  variable  treatment  is  valid. 
Additional  evidence  for  this  observation  is  discussed  in  Section  5.  In  stress 
relaxation  testa  (sponsored  by  LPC)  of  the  propellant  formulation  identical  to 
that  discussed  in  Section  3.3,  but  made  approximately  one  year  later  from 
different  lots  of  raw  materials,  the  viscoelastic  spectral  properties  and  shift 
factor  variation  with  temperature  were  observed  to  be  in  good  engineering 
agreement  with  the  data  presented  in  this  report.  This  type  of  agreement, 
while  not  as  yet  supported  statistically,  argues  effectively  for  the  basic  valid¬ 
ity  of  the  viscoelastic,  reduced  variable  approach  in  engineering  analysis  of 
propellant  grains. 
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Section  4 

GRAIN  THERMAL  STRESS  ANALYSIS 


4.  1  BACKGROUND 


Methods  are  well  established  for  analysis  of  the  quasl-statlc  boundary  value 
problem  for  homogeneous.  Isotropic,  linear  viscoelastic  materials  under 
surface  tractions,  body  forces  and  isothermal  temperature  conditions. 

Within  the  scope  of  small  deformations  and  time -Independent  boundaries,  the 
solutions  can  be  obtained  by  application  of  the  Laplace  transform,  which 
eliminates  time  dependence  and  reduces  the  viscoelastic  problem  to  an 
associated  elastic  one  (Lee,  1955). 

Introduction  of  either  a  transient  or  steady-state  thermal  field  produces 
difficulties  in  the  solutions  that  are  not  easily  surmounted.  The  material 
behavior  is  sensitive  to  temperature,  resulting  in  a  problem  where  the 
stresses  associated  with  the  thermally  induced  strains  are  sensitive  functions 
of  the  spatial  variables  and  time.  The  specific  nature  of  this  problem  was 
considered  in  subcontract  research  carried  out  by  Purdue  University  and 
methods  were  developed  for  obtaining  transient  thermal  stresses  in  hollow 
case-constrained  viscoelastic  cylinders.  The  details  of  the  work  are  re¬ 
ported  in  Volume  II,  Chapter  1  of  this  report. 

Two  methods  of  solution  developed,  and  the  one  best  suited  for  computer 
calculation  was  selected  for  numerical  solution  and  programed  for  the  IBM 
7090  computer.  This  solution  (description  begins  on  p.  35,  Volume  II)  con¬ 
siders  a  case -constrained,  circular  port  grain  under  plane  strain  end  con¬ 
straint.  The  propellant  was  assumed  to  be  incompressible  and  the  case 
elastically  rigid.  The  thermal  fields  were  assumed  to  be  axisymmetric  and 
externally  imposed.  Thermal  field  specification  and  propellant  viscoelastic 
properties  were  introduced  as  indicated  in  Sections  2  and  3  of  this  volume. 

s 

Various  types  of  thermal  inputs  specified  at  the  grain-case  interface  were 
examined  numerically.  The  possibility  of  stress  accumulation  was  investi¬ 
gated  for  various  step-cycle  thermal  Inputs  and  parametric  data  were  com¬ 
puted  showing  the  variation  of  peak  stress  with  grain  web  fraction,  grain 
diameter  and  temperature  for  step  surface  temperature  input. 

The  above  work  is  discussed  in  Section  4.  2  and  numerical  data  are  presented. 
Section  4.  3  discusses  the  possibility  of  estimation  of  peak  thermal  stresses 
using  simplified  limit  engineering  techniques.  Subcontract  effort  in  experi¬ 
mental  photothermovlscoelasticity  is  reviewed  in  Section  1.4.  Discussion  of 
the  thermal  stress  problem  and  the  engineering  significance  of  the  data  is 
presented  in  Section  4.  5. 


4-1 


LOCKHEKD  PROPULSION  COMPANY 


578-F 


4.  2  TRANSIENT  THERMAL  STRESS  CALCULATIONS 


The  transient  thermal  stress  calculations  were  based  on  the  following  set  of 
equations,  which  are  summarized  here  from  Sections  2  and  3  of  this  volume, 
and  from  Chapter  I  of  Volume  II.  The  symbols  used  are  defined  in  Table  4>1. 

Thermal  Field  Equations 

The  following  expressions  are  defined  (Ref.  Section  2.  2,  Vol.  I) 

O  ;  (1-i^  )  >  Ri  (t) 

Az  =  [  1  -  qz(T)  1  [1  -  ]*  +  T^Ct) 

where 

qi(T)  =  V  11-0299  T 

T  (T)  =  1  -  e'«**”*rt/-  -1) 
a  ta 


(4-1) 

(4-2) 

(4-3) 

(4-4) 


F  =  nondimensional  time  =  ri* 
t^  =  nondimensional  penetration  time  = 


Hence,  the  complete  step  solution  for  the  thermal  field  is 
A(r,r  )  =  Ti  Ai  (r,t)  1  (t^-T)  +  A2(r,T  )  1  (T- 1^) 


where 

1  (Fa  -T)  = 

1  (T  -FJ  = 


1.0  ;  t  <  ta 
O  ;  T  >  t^ 
O  ;  T<  ta 
1.0  ;  T>  ta. 


Ti  =  step  change  on  outer  case 

X  =  b/a 

Y  =  r/b 

F  =  K  t/b* 


(4-5) 
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A(r,  t) 


Ai(r.t) 

Ai(r,t) 

F(r.t) 

G(t)  =  E^(t) 
T(r.t) 

T 

a 

b 


Table  4-1 

LIST  OF  SYMBOLS 

Temperature  distribution  for  a  unit  step  change  on  the  grain 
outer  surface 

Function  defined  by  Equation  (4-1) 

Function  defined  by  Equation  (4-2) 

Function  defined  by  Equation  (4-10) 

Relaxation  modulus  of  propellant 

Temperature  change 

Temperature  change  on  outer  surface 

Inner  radius  of  grain 

Shift  factor 

1/a^ 

Outer  radius  of  grain 
Thermal  diffusivity 

Penetration  depth  of  heat  front,  Equation  (2-3) 

Function  defined  by  Equation  (2-16) 

Radial  coordinate 

Nondimensional  radial  coordinate  =  r/b 
Time 

Nondimensional  time  =  nt/b* 

Penetration  time  of  heat  front  to  r  =  a 
Function  defined  by  Equation  (4-8) 

Linear  coefficient  of  thermal  expansion 
Parameter  defined  by  Equation  (2-17) 

Hoop  strain 
Radial  strain 
b/a 

Function  defined  by  Equation  (4-11) 

Hoop  stress 
Radial  stress 

Specific  heat 
Thermal  conductivity 
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The  temperature  field  for  a  series  of  step  changes  on  the  exterior 
surface  is  given  by 

T(7.T)  =  r  T,n  A(7.t-tn)  1  (T-tJ  (4-6) 

where 


t^  =  1,  2,  3,  4  .  .  .  n  =  time  at  which  the  step  change  is  applied. 
Tin  =  step  change  in  the  surface  temperature  at  t  =  tn- 
Viscoelastic  Analysis 


The  equations  for  the  viscoelastic  analysis  are 
(Ref.  Chap.  1,  Vol.  II,  Section  3.3,  Vol.  I) 


€g(r,  t)  =  (r  ,  t) 

(4-7) 

_  ,  F  _ 

'l'(r.t)  =  -y4jfpT(p,t)dp 

(4-8) 

T 

Oj.  (r,  T  )  =  J  F  (p  ,  t)d  p 

(4-9) 

whe  re  _ 

f(7,T)  =  ^  |g[C(T,7)  -  c(t,7)]|-^  [2¥(7,t).t(7,  T)]dT 
C  =  Y  J  ^  r  T  (7,  T  )  1  dx 

G{t)  =  Ej.{t)  =  330t"®'*‘’ 


(4-10) 

(4-11) 

(4-12) 


The  function  i  (T)  is  the  reciprocal  of  a|^  (Section  3.3).  The  hoop 
stress  is  given  by 


^  =  F(7,  T)  (4-13) 

Sr 

7-^^  +  or  =  0p(7,  T  )  (4-14) 

Sr  ^ 

and  the  radial  strain  specified  by 
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Hoop  Stress  at  Port  versus  Time  for  Step,  Step  with  a  Jog,  and  Cyclic 
Thermal  Inputs 
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on  the  peak  stresses  obtained  as  a  result  of  step  environment  changes  in  a 
one -foot  diameter  motor.  The  calculations  considered  the  following  para¬ 
metric  variations; 

•  Thermal  step  inputs  of  -40®F.  -100®F.  -130®F,  and  -160®F 
from  a  reference  temperature  of  70®F. 

•  Web  fractions  of  9,  50,  66.7,  75  and  80  percent. 

The  resulting  peak  stresses  are  shown  in  Figures  4-5  through  4-8. 

Peak  stresses  for  various  motor  diameters  can  be  determined  from  the  data 
in  Figures  4-5  through  4-8  by  the  following  scaling  law: 

•  Obtain  the  stresses  Oq  for  a  given  X  (grain  outer  diameter  to 
inner  diameter  ratio)  and  a  particular  grain  outer  diameter,  Zbo. 

•  For  a  new  grain  outer  diameter  =  Zb^,  the  new  stress  value  Oy 
is  obtained  as 

o  1  =  }  0  0 

This  scaling  law  is  applicable  only  for  the  particular  viscoelastic  properties 
of  the  Polycarbutene  R  propellant  studied  in  this  program.  In  general,  for 
any  power  law  approximation  to  the  relaxation  modulus  of  the  form 


Er  (t)  =  K  t'*" 


the  exponent  of  the  approximation,  m,  can  be  used  to  obtain  the  scaling 
law  as 


The  ratio  ai/oo  for  various  motor  diameter  ratios  is  shown  in  Figure  4-9. 
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GRAIN  O.D./I.D.  RATIO,  \ 


Figure  4-7  Case-Grain  Interface  Pressure  versus  Temperature 
Step  and  Grain  O.  D.  to  I.  D.  Ratio 
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Figure  4-8  Case -Grain  Interface  Pressure  versus  Temperature 
Step  and  Grain  O.  D.  to  I.  D.  Ratio 
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4.  3  ELASTIC  AND  VISCOELASTIC  LIMIT  ANALYSIS 


In  many  instances  it  is  useful  to  calculate  limits  for  stresses  in  a  body  on 
the  basis  of  assumptions  that  enable  simplification  in  computation,  but  which 
lead  to  reliable  values  for  limits  within  which  the  rigorously  derived  answer 
must  lie.  \In  the  thermal  transient  problem  the  following  two  approaches  are 
obvious. 

•  Calculation  of  the  limit  stresses  associated  with  the  glassy 
and  rubbery  propellant  moduli. 

•  Calculation  of  the  stress  assuming  a  step  function  body 
change  in  the  grain  temperature,  using  the  propellant 
properties  at  the  new  temperature  and  solving  for  the 
time  at  which  thermal  equilibrium  is  attained. 

Noting  that  the  glassy  and  rubbery  equilibrium  stresses  will  be  different  by 
approximately  three  orders  of  magnitude,  the  first  approach  can  be  ruled 
out  without  detailed  consideration.  Equilibrium  stresses  for  the  rubbery 
condition  are  obtained  only  after  long  times  at  moderate  temperatures  and 
at  high  temperatures  after  short  times.  The  thermal  stress  calculations 
indicate  values  for  these  stresses  of  the  order  of  1  psi  for  temperatures 
above  70*F  for  Polycarbutene  R  propellant  for  the  more  severe  cases  (small 
motor  diameter,  high  web  fraction).  Equilibrium  stresses  associated  with 
glass -like  behavior  are  physically  unattainable  even  by  shock  cooling  of  a 
real  grain  because  of  finite  thermal  equilibrium  times  and  the  virtual  im¬ 
possibility  of  cooling  a  practical  grain  below  its  glass  transition  temperature 
without  failure  by  web  cracking. 

The  second  approach  is  more  promising,  but  is  arbitrary  in  that  a  thermal 
equilibrium  time  must  be  assumed.  The  thermal  equilibrium  time  is  mathe¬ 
matically  infinite  and  therefore  does  not  permit  selection  of  a  non-arbitrary 
time  for  solution.  Taking  the  time  to  be  that  at  which  the  grain  reaches 
98  percent  of  thermal  equilibrium  (about  the  time  that  the  peak  stresses  are 
obtained  in  the  transient  step  temperature  solution)  leads  to  a  calculated 
core -grain  interface  stress  of  15  psi  for  a  1-foot  diameter,  50  percent  web 
fraction  grain.  This  compares  to  the  computer  value  of  about  9  psi  (Figure  4-2). 

While  in  this  instance  the  limit  calculation  yields  an  estimation  for  the  stress 
that  is  in  the  "ballpark",  the  size  of  the  ballpark  can  be  expected  to  change 
an  analytically  unpredictable  amount  with  temperature,  grain  diameter,  and 
web  fraction.  The  calculation  yields  stress  values  that  are  probably  better 
than  an  uninformed  guess,  but  of  little  real  engineering  use. 
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4.4  PHOTOTHERMOVISCOELASTICITY 


The  use  of  photothermoviscoelastic  techniques  as  an  experimental  adjunct  to 
transient  thermal  stress  analysis  was  studied  in  this  program,  largely  in 
subcontract  activities  at  the  University  of  Washington  and  at  New  York  Uni¬ 
versity.  The  purpose  of  the  New  York  University  program  was  to  obtain, 
variously,  photoviscoelastic  data  for  transient  thermal  fields.  The  Univer¬ 
sity  of  Washington  effort,  much  smaller  in  scope,  was  directed  at  theoreti¬ 
cal  inquiry  in  this  area  of  interest. 

The  New  York  University  effort,  under  the  personal  direction  of  H.  Becker, 
proceeded  in  a  photoviscoelastic  calibration  of  a  Hysol  material  at  tempera¬ 
tures  near  the  glass  transition  temperature  of  the  material.  Photographs 
of  fringe  patterns  in  simulated  plane  stress  rocket  grain  configurations 
made  from  the  Hysol  were  also  obtained.  However,  the  work  did  not  yield 
information  of  specific  technical  worth  to  the  program, 
largely  because  of  insufficient  development  of  the  theoretical  aspects  of  the 
problem*. 

Redirection  of  the  photothermoviscoelastic  research  was  made  to  permit 
concentration  on  the  theoretical  problems,  with  the  redirected  work  being 
carried  out  at  the  University  of  Washington  under  the  direction  of  Prof. 

R.  J.H.  Bollard  and  with  the  assistance  of  Dr.  Ellis  Dill.  A  report  of  this 
investigation  is  presented  in  Section  4.  4.  1 


‘The  final  report  of  N.Y.U.  effort  was  submitted  on  12  Nov  1962  with 
LPC  Monthly  Progress  Report  No.  11. 
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4.  4.  1  University  of  Washington  Subcontract  Work* 


Introduction.  A  theoretical  basis  was  developed  for  the  experimental  deter- 
mlnation  of  the  stress  state  of  a  viscoelastic  material  in  transient  plane 
stress  due  to  temperature  changes  and  external  loads. 

The  following  summarizes  the  theoretical  and  experimental  results  of  this 
work. 


Basic  Theory.  Experiments  on  the  birefringence  of  high  polymers  indicate 
that  the  birefringence  depends  on  the  time  history  of  the  mechanical  state  in 
the  same  way  as  constitutive  relations  for  linear  viscoelastic  materials. 

This  suggests  the  hypothesis  that  the  dielectric  tensor  to  the  one-half  power 
depends  on  the  strain  tensor  by  a  Duhamel  type  or  hereditary  integral,  of 
the  same  form  as  for  the  dependence  of  the  stress  tensor  on  strain  history, 
but  with  a  different  kernel  which  may  be  called  the  optical -relaxation function. 
For  isothermal  deformations  this  leads  to  a  relation  between  stress  and  fringe 
order  n  of  the  form 

^x  ■ '^y  ■  K  J I  ^(  ?)  cos  2a  {? )  I  d§  , 

txy  =  ^  J  K(t  -  I  il(§)  sin  2a  (§)|d?  ; 

where  h  is  the  model  thickness  and  a(t)  is  the  angle  between  the  x-axis  and 
the  axis  of  polarization. 

If  it  is  found  that  the  temperature  dependence  of  the  optical -relaxation 
function  can  be  described  by  a  temperature  shift-factor  which  is  identical 


*This  section  was  prepared  by  Professor  R.  J. H.  Bollard  and  Associate 
Professor  E.H.Dill,  University  of  Washington. 
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to  that  for  the  mechanical  response,  then  the  thermal  stresses  depend  on 
the  birefringence  in  a  similar  way*: 

'f 

Ox  -  ^  J*  K  (t-  ^  I  Ti  (?)  COB  2a  {?)|d?  , 


"^y  =  ^  J*  k(t-  ?)^  I  ri(?)  sin  2a{?)|d?  ; 

O  "  . 

where  t  (t)  is  the  reduced  time, 
t 

T(t)  =  J  $(T(§)d?  ; 
o 

$(T)  is  the  temperature  shift  factor; 

)*  ; 


and 


K(t)  =  K($(T)t) 

is  the  optical  relaxation  function  at  temperature  T. 


Experimental  Procedure.  According  to  the  theory,  it  is  necessary  to  have 
a  complete  record  of  fringe  order  n(t),  the  angle  of  polarization  a(t),  and 
the  temperature  T(t)  at  a  point  in  order  to  determine  the  difference  in  prin¬ 
cipal  stresses  and  the  axis  of  principal  stress  at  that  point. 

The  fringe  order  is  found  by  observation  of  the  isochromatic  lines.  A  number 
of  photographs  of  the  isochromatic  pattern  at  a  sequence  of  time  s 
ti,  tz . tj^  will  allow  the  determination  of  n(t)  at  each  point. 

In  order  to  find  the  axis  of  polarization  at  an  arbitrary  point  it  is  necessary 
to  make  simultaneous  observations  of  the  isoclinic  pattern  for  various  angles 
of  polarization  of  the  incident  light.  In  an  ordinary  photoelastic  polariscope 
this  is  accomplished  by  manually  rotating  the  polarizing  and  analysing 
element.  However,  in  the  transient  viscoelastic  problems  which  may  be  of 
interest  the  isoclinic  field  may  vary  rapidly,  and  it  will  be  necessary  to 
make  perhaps  ten  photographs  of  the  isoclinic  field  for  ten  different  orienta¬ 
tions  of  the  polarizer  and  analyser.  The  need  for  this  "instantaneous"  sweep 
of  the  isoclinic  field  precludes  the  use  of  manually  rotated  elements  in  the 
bench.  Mr.  Michael  Fourney  has  suggested  that  the  rapid  rotation  of  the 
angle  of  polarization  of  the  incident  light  can  be  accomplished  by  means  of 
‘Appendix  B 
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the  Faraday  effect.  This  electronic  means  will  allow  the  recording  of  ten 
photographs  of  the  isoclinics  for  each  photograph  of  the  isochromatic  pattern 
at  practically  the  same  time  instant  ti. 

The  photographic  record  can  then  be  analyzed  to  determine  n(t)  and  a(t)  for 
any  point. 

The  complete  experimental  apparatus  is  diagramed  in  Figure  4-10.  The 
details  of  operation  for  a  bench  of  this  type  were  described  in  the  earlier 
report. 


Program  for  Experiment.  Phase  one  of  the  experimental  program  should  be 
to  determine  whether  there  exists  model  materials  which  have  the  desired 
mechanical  properties  and  for  which  a  linear  integral  relation  between  strain 
tensor  and  dielectric  tensor  exists.  Some  initial  tests  have  been  carried  out 
in  this  program  to  determine  whether  further  studies  are  justified. 

A  tension  test  creep  specimen  was  placed  in  a  standard  plane  polariscope. 

The  change  in  length  was  recorded  automatically  as  a  function  of  time  by 
means  of  an  LVDT  displacement  gage.  The  fringe  order  variation  with  time 
was  determined  from  continuous  recording  of  the  intensity  of  light  transmitted 
by  the  model  with  a  photocell.  The  apparatus  is  shown  in  Figure  4-11. 

The  most  promising  material  found  was  a  polyester -styrene  resin  with  trade 
name  Polylite  8151  supplied  by  Reichold  Chemical  Company  to  which  1^  MEK 
peroxide  hardener  was  added  and  the  batch  poured  in  the  mold  and  cured  at 
room  temperature.  The  results  of  the  test  are  purely  qualitative  and  are 
shown  in  Figures  4-12  and  4-13. 

It  can  be  seen  that  the  optical  creep  is  a  curve  of  the  same  form  as  the 
mechanical  creep,  but  not  exactly  the  same.  Thus,  the  optical  creep  function 
can  not  be  a  simple  scalar  factor  times  the  mechanical  creep  compliance. 

A  more  elaborate  experimental  calibration  apparatus  has  been  constructed 
which  will  allow  an  accurate  determination  of  the  limits  of  applicability  of 
the  linear  integral  relation  at  various  temperatures,  and  to  investigate  the 
existence  of  an  optical  temperature  shift  factor  for  the  Polylite  and  other 
materials.  This  apparatus  is  diagramed  in  Figure  4-14.  This  bench  has 
recently  been  completed  and  is  shown  in  Figure  4-15. 

Phase  two  of  the  experimental  program  will  be  to  develop  methods  of  practical 
stress  analysis.  In  conjunction  with  this  phase  the  experimental  apparatus  of 
Figure  4-10  will  be  assembled.*  It  is  hoped  to  have  an  operating  bench  for 
quantitative  investigations  by  midsummer  of  1963. 


*Partially  supported  by  a  grant  from  the  National  Aeronautics  and 
Space  Administration. 
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Figure  4-10  Plane  Polariscope  for  Transient  Loadings 
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Figure  4-11  Mechanical  and  Optical  Creep  Test 
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Figure  4-14  Diffused  Light  Polariscope  for  Material  Characterizations 
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Figure  4-15  Material  Calibration  Apparatus 
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RecommendationB .  Steps  should  be  taken  to  establish  at  least  two  different 
centers  to  undertake  experimental  characterization  of  materials  which  may 
be  used  to  mechanically  simulate  solid  propellant  materials  and  which  ex¬ 
hibit  the  birefringent  effects.  It  appears  that,  with  the  numerous  resins 
and  plasticizers  available,  mechanical  properties  of  nearly  arbitrary  nature 
can  be  achieved.  It  is  recommended  that  some  standard  material  and  curing 
procedure  be  agreed  upon  by  the  various  centers.  The  experimental  program 
should  then  proceed  in  the  indicated  manner  --  first,  complete  characteriza¬ 
tion  of  the  standard  material;  second,  development  of  a  routine  stress  analysis 
technique. 

Theoretical  work  toward  establishing  a  basis  for  applying  photoanalysis  where 
nonlinear  viscoelastic  behavior  is  present  should  be  started  at  this  time.  It 
is  felt  that  photoanalysis  may  prove  to  be  an  effective  tool  for  such  problems. 
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4. 5  DISCUSSION 

4.  5.  1  Transient  Thermal  Stress  Computation 


The  data  presented  in  Section  4.  2  for  transient  thermal  stresses  in  a  circular 
port,  plane  strain,  rigid  case  propellant  grain  are  the  only  such  data  known  to 
be  available  as  of  this  date.  The  engineering  validity  of  the  data  for  thermal 
stresses  is  subject  to  unknown  but  presumably  small  error  limits.  Following 
is  an  examination  of  the  transient  analysis  data,  considering  the  accuracy  of 
the  various  input  parameters  to  the  problem,  and  a  review  of  the  engineering 
applicability  of  the  analysis. 

Thermal  field  input  information  as  entered  in  the  computer  analysis  explicitly 
determines  the  grain  strains  and  the  spatial  dependence  of  the  grain  properties 
Hence,  the  associated  stress  field  is  that  required  by  the  thermal  field  and 
its  interaction  with  the  grain.  As  specified  by  the  field  expressions  discussed 
in  Section  2,  Volume  I,  thermal  distributions  have  been  verified  to  be  accurate 
for  specified  temperature  boundary  conditions  both  in  motors  and  in  uncured 
propellant  grains.  In  both  the  particular  case  of  a  grain  surface  (or  motor 
case)  temperature  specification  and  in  the  general  case  of  a  specification  of 
a  motor  environmental  temperature,  the  thermal  field  solutions  for  the  grains 
converge  in  accuracy  as  thermal  equilibrium  is  approached.  The  sensitivity 
of  the  stresses  is  highest  to  thermal  field  variation  near  the  thermal  equilib- 
ruim  condition  for  a  negative  temperature  change.  Therefore,  the  thermal 
field  predictions  are  estimated  to  be  of  excellent  engineering  accuracy  for 
a  given  thermal  input  in  which  the  case  temperature  is  decreased  in  some 
time -dependent  manner,  arriving  at  a  specific  constant  level.  The  above 
comments  are  predicated  upon  axisymmetric  thermal  field  conditions. 


The  propellant  viscoelastic  data  used  in  the  thermal  stress  problem  are 
probably  the  likeliest  source  of  error.  This  is  partially  because  of  an 
inherent  difficulty  in  obtaining  broad  spectrum  physical  data  for  a  propellant 
and  certain  pragmatic  objections  to  the  reduced  variable  approached  used  in 
construction  and  application  of  the  master  viscoelastic  representation.  The 
master  or  reduced  relaxation  specification  (Section  3.2,  Volume  I)  is  con¬ 
structed  from  isothermal  tests  across  discrete  time  intervals  for  various 
temperatures.  Hence,  the  master  relaxation  spectrum,  which  extends 
across  approximately  20  decades  of  time  at  a  single  temperature,  is  ass¬ 
embled  from  segments  of  isothermal  data  each  of  which  is  rigorously  valid 
only  across  the  2  to  5  decades  of  time  of  actual  observation.  The  associated 
time -temperature  shift  relation,  aq-,  is  developed  as  a  product  relation  in 
the  reduction  procedure,  a  technique  inherently  susceptible  to  cumulative 
inaccuracy.  The  accuracy  of  the  shift  factor,  a-p,  is  of  paramount  impor¬ 
tance  in  the  thermal  stress  calculations  since  the  sensitivity  of  the  stresses 
to  variation  of  a-p  in  the  solution  is  higher  than  for  any  other  single  parameter 
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Additionally,  the  fundamental  question  arises  as  to  the  validity  of  equating, 
for  example,  a  160®F  temperature  change  to  a  10-decade  time  scale  shift. 
Theory  in  this  connection  is  not  particularly  well  developed;  however, 
experimental  observations  do  lend  credence  to  the  results  obtained  in  the 
reduced  variable  approach  (Section  5,4,  Volume  I). 

The  preceding  commentary  on  the  accuracy  of  the  propellant  viscoelastic 
properties  presents  a  pessiinistic  view  of  the  problem,  but  probably  an 
unrealistically  pessimistic  view.  It  was  noted  previously  that  a  replicate 
viscoelastic  characterization  of  the  Polycarbutene  R  propellant  conducted 
on  a  batch  cast  from  different  lots  of  raw  materials  gave  reduced  properties 
essentially  identical  to  the  data  used  in  this  program  (Section  3.5,  Volume  I). 
The  observation  shows  that  the  viscoelastic  data  are  generally  reproducible, 
a  vital  intrinsic  characteristic  for  engineering  analysis.  In  addition,  the 
nature  of  the  stress  analysis  problem  is  such  that  errors  arising  from  pro¬ 
pellant  property  error  would  tend  to  increase  with  the  magnitude  of  the 
temperature  excursion  from  the  reference  state.  Restricting  the  tempera¬ 
ture  range  of  interest  to  -65®F  and  above,  the  errors  appear  tolerable  from 
the  point  of  view  that  the  computed  stress  fields  are  consistent  with  observed 
grain  failure  behavior  (Section  5.4,  Volume  l). 

In  summary,  while  determination  of  the  absolute  accuracy  of  propellant 
reduced  variable  data  for  stress  analysis  purposes  requires  additional 
sophisticated  experimental  inquiry,  acceptable  engineering  accuracy  in 
propellant  broad  spectrum  properties  characterization  has  been  attained, 
as  judged  by  available  criteria. 

The  asumptions  made  in  the  transient  thermal  stress  analysis  are  reasonable 
engineering  idealizations  of  practical  rocket  motor  grains.  Stipulation  of 
propellant  rubberlike  incompressibility  has  been  shown  to  be  an  accurate 
engineering  condition  in  independent  investigations  (Jones,  1963).  The 
assumption  enables  certain  analytic  simplifications  in  the  stress  analysis 
solution,  among  which  is  the  simplification  of  the  basic  stress -strain  law 
for  rigorous  interpretation  in  a  state  of  uniaxial  tension,  the  physical  state 
in  which  the  propellant  properties  were  measured.  The  assumption  of  plane 
strain  end  conditions  for  the  grain  tends  in  the  conservative  direction  for 
most  moderate  length,  moderate  web  fraction  motors.  In  general,  a  plane 
strain  analysis  that  predicts  undesirably  high  stress  or  strain  levels  in  a 
propellant  grain  signals  a  need  for  extreme  care  in  the  design  of  a  finite 
length  grain  for  the  same  diameter  and  web  fraction.  The  assumption  of 
case  rigidity  is  consistent  with  the  observed  behavior  of  steel  case  rocket 
motors  under  low  states  of  grain  stress,  such  as  are  observed  under  trans¬ 
ient  thermal  states. 

The  basic  analysis  as  developed  by  Purdue  University  was  shown  to  have  an 
analytic  proof  of  convergence.  The  accuracy  of  the  calculations  inherent  to 
the  mathematic  framework  was  shown,  hence,  to  be  restricted  only  by  the 
extent  to  which  temporal  iterations  and  calculation  of  series  expansion 
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terms  was  done  (Volume  II,  Appendix  III).  Both  of  these  latter  processes 
are  limited  by  the  computer  storage  capacity  and  are  practical  rather  than 
theoretical  problems.  In  practice,  the  observed  reproducibility  of  the 
computer-determined  stresses  was  approximately  ±  2  percent,  the  limita¬ 
tion  on  accuracy  arising  from  the  temporal  iteration  flexibility  of  the  program. 

In  summary,  the  stresses  obtained  in  the  calculation  are  believed  to  be  of 
realistic  engineering  value  for  the  imposed  thermal  boundary  conditions. 


4.  5.  2  Thermal  Stress  Accumulations 


A  primary  objective  in  this  program  was  elucidation  of  the  possible  phenomenon 
of  thermal  stress  accumulation.  Theoretical  study  of  thermal  stresses  in  a 
viscoelastic  slab  arising  from  step  cyclic  temperature  boundary  conditions 
showed,  for  uncomplicated  constraints,  that  thermal  stress  accumulation  was 
conceptually  possible.  It  was  not  possible  during  the  program  to  develop 
analytic  proof  or  disproof  for  the  phenomenon  for  the  case -bonded  circular 
port  grain  subjected  to  external  thermal  inputs.  Thus,  the  numerical  data 
for  various  types  of  step-cyclic  external  thermal  inputs  provide  the  only  in¬ 
sight  into  the  problem  of  stress  accumulation  in  rocket  grains.  The  data  show 
no  evidence  for  stress  accumulation  in  a  linear  viscoelastic  propellant  grain, 
as  considered  in  Section  4.  2,  Volume  I,  for  widely  varying  thermal  field  input 
conditions.  The  conclusion  reached  was  that  stress  accumulation  in  linearly 
viscoelastic  propellant  grains  under  realistic  field  conditions  is  unlikely  in 
any  significant  degree,  even  though  the  conceptual  possibility  for  accumula¬ 
tion  represents  an  as  yet  unresolved  problem 


4.  5.  3  Peak  Thermal  Stresses  for  Shock  Cooling 


The  parametric  data  for  peak  thermal  stresses  in  grains  subjected  to  various 
negative  thermal  step  inputs  (Figures  4-5,  6,  7,  and  8)  was  tabulated  to  pro¬ 
vide  a  basis  for  engineering  limit  analysis.  Specifically,  the  stress  levels 
shown  correlate  to  the  maximum  stresses  obtained  following  instantaneous 
thermal  coupling  of  the  rocket  motor  case  to  an  infinite  negative  heat  source 
(e.g.  ,  immersion  of  the  motor  in  a  fluid  at  the  reduced  temperature,  for 
engineering  purposes).  Higher  stress  peaks  are  not  possible  for  any  physically 
realistic  process  by  an  external  cooling  process  without  stress  accumulation. 
The  peak  stresses  shown  are  conservatively  high  for  slower  rates  of  cooling 
and  therefore  represent  useful  criteria  for  engineering  grain  analysis.  Con¬ 
siderations  regarding  grain  failure  criteria  and  structural  integrity  analysis 
are  discussed  in  Section  5,  Volume  I,  and  in  Volume  III  of  this  report. 
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4.5.4  Photothermoviscoelasticity 


During  this  program,  the  state  of  the  art  in  use  of  photoviscoelastic  investi¬ 
gation  of  thermal  stresses  regressed  from  a  condition  of  engineering  appli¬ 
cability  (as  originally  proposed  by  New  York  University)  to  a  state  of  con¬ 
ceptual  theoretical  study.  It  was  found  upon  rigorous  investigation  of  the 
problem  that  observed  fringe  patterns  in  transiently  heated  viscoelastic 
grain  models  could  not  be  interpreted  for  quantitative  stress  values  because 
of  an  insufficient  theoretical  and  experimental  basis  for  the  interpretation. 
However,  theoretical  and  experimental  study  at  the  University  of  Washington 
resulted  in  significant  progress  in  this  area.  A  theoretical  foundation  was 
developed  for  interpretation  of  stresses  in  birefringent,  viscoelastic  models, 
and  limited  experimental  progress  was  made  in  verification  of  the  theory. 

Continued  effort  in  development  of  experimental  stress  analysis  techniques 
for  the  transient  thermal  viscoelastic  problem  is  considered  important.  At 
the  present  time,  the  possibility  of  analytic  solution  for  transient  thermal 
stresses  in  a  star-perforated  grain,  or  in  a  grain  asymmetrically  heated, 
is  considered  unlikely  by  presently  known  techniques.  While  a  problem  of 
considerable  difficulty,  the  development  of  photothermoviscoelastic  experi¬ 
mental  techniques  may  be  tractable.  If  so,  engineering  treatment  of  struct¬ 
ural  problems  particularly  pertinent  to  high-performance  airborne  missiles 
may  prove  feasible. 
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Section  5 

GRAIN  FAILURE  ANALYSIS 


5.  1  BACKGROUND 


A  conoplete  grain  structural  integrity  analysis  requires  the  existence  of  a 
suitably  accurate  method  for  determining  grain  stresses  and  strains  and  a 
criterion  for  evaluation  of  the  possibility  of  failure  of  the  grain.  Grain  stress 
analysis  for  isothermal  grain  conditions  of  both  general  and  specific  interest 
presently  is  in  an  advanced  state  for  linear  materials  under  small  strain  con¬ 
ditions  (Fitzgerald,  1962).  The  state  of  the  art  in  engineering  analysis  for 
transient  thermal  states  is  updated  by  this  report.  However,  propellant  grain 
failure  criteria  for  either  isothermal  or  transient  thermal  states  are  some¬ 
what  more  tenuous  on  both  theoretical  and  experimental  bases. 

This  section  presents  the  results  of  experimental  grain  failure  investigations 
for  transient  and  cyclic  thermal  conditions  as  investigated  in  scaled  or  analog 
rocket  motors. 
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5.  2  MOTOR  ANALOG  EXPERIMENTS 


A  set  of  4-inch  diameter  by  14-inch  long,  circular  port  analog  rocket  motors 
(Figure  5-1)  of  varying  web  fraction  were  cast  from  a  single  batch  of  Poly- 
carbutene-R  propellant.  The  dimensional  data  for  the  motors  are  as  follows: 

Length  Grain  Outer  Diameter  Web  Fraction 
Motor  No.  T  Diameter  Grain  Inner  Diameter  (percent) 


lA,  IB 

2.5 

4/1.5 

63.5 

2A,  2B 

2.5 

4/1.2 

70 

3A,  3B 

2.5 

4/1.0 

75 

2.5 

4/0.866 

78 

The  motors  were  subjected  to  various  externally  imposed  thermal  conditions 
with  the  motors  sealed,  internally  desiccated,  and  with  the  ends  of  the  grain 
insulated.  All  motors  were  stored  in  a  70®  ±3°F,  desiccated  environment  for 
3/2  months  prior  to  initiation  of  the  test  series.  The  motor  case  and  port 
surface  temperatures  were  monitored  with  bonded,  standardized  thermo¬ 
couples  in  all  tests.  Port  diameter  measurements  in  the  motors  were  made 
at  the  midpoints  of  the  grains  using  a  Mueller  gage  inserted  through  a  tem¬ 
porarily  unsealed  end  of  the  motor  (Figure  5-2).  In  all  cases  where  strain 
measurements  were  made,  the  operation  was  carried  out  quickly  enough  to 
confine  the  grain  thermal  excursion  during  measurement  to  less  than  1®F. 
Grain  port  measurements  were  reproducible  to  within  ±0.002  inch. 


5.  2.  1  Rapid  Monotonic  Cooling  Test 


Motors  No.  1,  2,  3  and  4  were  subjected  to  rapid  cooling  from  equilibrium  at 
70®F  to  equilibrium  at  -75®F,  Figures  5-3,  5-4,  5-5,  and  5-6  show  the  case 
and  port  temperatures,  and  the  measured  port  hoop  strain,  calculated  as  the 
incremental  change  in  diameter  divided  by  the  original  casting  mandrel 
diameter.  Figure  5-5  shows  the  predicted  grain  port  temperature  versus 
time  based  on  the  techniques  discussed  in  Section  2,  Volume  1,  using  the 
observed  case  temperature  as  the  pertinent  boundary  condition. 

Based  on  the  cumulative  strain  at  maximum  stress  technique  (Jones,  1963), 
grain  failure  predictions  were  made  taking  the  temperature  and  strain  rate 
conditions  of  pertinence  to  be  those  measured  at  the  grain  ports. 
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Figure  5-1  5 -inch  Diameter  by  14 -inch  Long  Circular  Port  Analog  Motor 
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LENGTH  TO  DIAMETER  RATIO  14/4 
GRAIN  O.  D./GRAIN  I.  D.  =  4/1.5 
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The  technique  is  as  follows.  For  a  condition  of  continuous  slow  cooling,  the 
port  hoop  strain  increments  for  small  temperature  increments  are  calculated. 
The  strain  occurring  in  each  increment  of  time  and  temperature  is  related  to 
the  propellant  physical  data  for  the  corresponding  strain  rate  and  temperature 
The  fractional  consumption  of  strain  capability  used  in  traversing  each  tem¬ 
perature  increment  is  svimmed.  The  temperature  at  which  the  sum  of  the 
fractions  equals  unity  is  the  predicted  low  temperature  limit  for  the  grain.  A 
description  of  the  analytical  and  interpretive  methods  is  presented  in  Volume 
III  of  this  report.  The  (isothermal)  strain  versus  temperature  profiles  for  the 
motors  and  the  predicted  failure  conditions  are  shown  in  Figure  5-7.  With  the 
exception  of  the  highest  web  fraction  analog  (No.  4),  the  failure  temperatures 
were  predicted  quite  accurately.  Failure  of  the  No.  4  motor  was  anticipated 
somewhat  below  the  isothermally  predicted  temperature  because  of  the  tran¬ 
sient  conditions.  The  tests  demonstrated  excellent  accuracy  in  prediction  of 
failure  temperatures  for  the  motors  of  further  interest,  the  replicates  of 
motors  No.  1,  2,  and  3. 


5.  2.  2  Thermal  Shock  Experiments 


Motors  No.  3A  and  3B  were  subjected  to  rapid  cooling  to  -30®F  and,  after 
obtaining  thermal  equilibrium,  immersed  in  a  circulating  water  bath  main¬ 
tained  at  193°F.  The  pertinent  thermal  field  data  (Figure  5-8)  show  excellent 
agreement  between  measured  and  predicted  port  temperature  versus  time 
variation.  The  predicted  port  temperature  variation  was  based  on  the  energy 
methods  described  in  Section  2,  Volume  I. 

With  reference  to  Figure  5-7,  at  -30®F  the  motors  were  within  15°F  (or  within 
3  percent  strain)  of  predicted  failure.  Limit  calculations  assuming  a  160®F 
temperature  rise  of  the  case  in  two  minutes  predicted  that  the  additional  port 
strain  increment  would  accumulate  the  additional  strain  to  the  failure  fraction 
required  to  cause  grain  cracking.  Alternately,  the  grain  stresses  also  were 
predicted  to  rise  to  values  near  those  associated  with  either  case  unbonding 
or  port  cracking  conditions.  The  motors  failed  by  extensive  port  cracking. 


5.  2.  3  Cyclic  Experiments 


Motors  No.  lA,  IB,  2A  and  2B  were  subjected  to  repetitive  cycling  between 
-50°F  and  160°F,  with  each  complete  cycle  taking  place  in  a  24-hour  period. 
(Figure  5-9).  With  reference  to  Figure  5-7,  at  -50®F  the  No.  2  motors  were 
within  10°  to  15°F  (or  3^2  percent  strain)  of  predicted  failure  and  the  No.  1 
motors  within  25°  to  30°F  (or  9  percent  strain)  of  failure.  One  of  the  No.  2 
motors  failed  after  one  cycle,  and  one  failed  after  eight  cycles.  The  No.  1 
motors  survived  22  cycles  without  sign  of  failure  and  without  change  in  the 
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Figure  5-7 


Analog  Motor  Experiments  Shock  Cooling 
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Figure  5-8  Shock'Heating  Experiment,  4  by  14-inch  Analog  Motors 
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Figure  5-9  Cyclic  Conditioning  Test,  4  by  14-inch  Analog  Motor 
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port  diameter  as  measured  at  thermal  equilibrium  at  70®F  during  each  cycle. 
Following  cycling,  the  No.  1  motors  were  cooled  rapidly  from  70®F  to  failure, 
which  was  observed  at  -80®  to  -85®F.  Hence,  it  was  observed  that  an  accumu¬ 
lation  process  leading  to  failure  took  place  in  the  No.  2  motors.  Conversely, 
the  No.  1  motors  exhibited  no  accumulation  of  failure  during  cycling  as  demon¬ 
strated  both  by  the  dimensional  stability  of  the  grains  and  the  final  failure 
temperatures  of  the  motors,  which  agreed  with  those  for  uncycled  motors 
(  reference  Figure  5-7). 

In  review  of  the  propellant  constant  strain  rate  uniaxial  data,  it  was  observed 
that  in  the  temperature  range  0®  to  -100®F,  apparent  nonlinearity  in  the  visco¬ 
elastic  behavior  was  observed  at  strain  levels  above  0.6  to  0.8  the  measured 
strain  at  maximum  stress  value,  cm  .as  reduced  for  time  and  temperature 
effects  .  By  comparison,  the  No.  1  motors  on  each  cycle  were  strained  to 
approximately  0.5  em  and  the  No.  2  motors  to  approximately  0.8  e  m,  at  -50®F. 

The  observations  are  consistent  with  the  hypothesis  that  the  strain  limit  defined 
by  linear  viscoelastic  behavior  of  the  propellant  corresponds  to  a  thermal 
cycling  fatigue  failure  limit.  While  not  statistically  imposing,  the  data  are 
vitally  interesting  and  indicate  an  important  area  for  further  experimental 
effort. 
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5.  3  DISCUSSION 


The  strains  in  a  transiently  cooled  grain  are  produced  only  by  the  thermal 
field  and  not  influenced  to  a  significant  degree  by  the  stresses.  Time  and  the 
temperature  integration  of  the  strains  and  strain  criteria  therefore  is  rela¬ 
tively  straightforward,  accepting  the  reduced  variable  treatment  of  the  uni¬ 
axial  failure  strain  data.  The  data  show  evidence  for  acceptance  of  the  lab¬ 
oratory  uniaxial  strain  at  maximum  stress  as  a  grain  failure  criterion  for 
monotonic  cooling  of  a  motor  to  failure.  Additional  support  for  this  criterion 
has  been  discussed  by  Jones,  1963.  With  some  trepidation,  the  limited 
experimental  motor  cycling  data  obtained  in  the  program  indicate  that  three - 
fourths  the  laboratory  strain  at  maximum  stress  is  a  reasonable  thermal 
cyclic  failure  criterion.  The  reduced  variable  approach  to  comparison  of  the 
strains  is,  of  course,  necessary. 

However,  evaluation  of  the  motor  analog  data  in  terms  of  possible  stress- 
determined  failure  criteria  was  not  possible.  Such  investigation  would  have 
required  theoretical  and  experimental  effort  outside  the  scope  of  the  program. 
Complex  by  comparison  to  evaluation  of  the  analog  data  by  strain  criteria, 
the  nature  of  the  problem  i  s  as  follows,  T  r  a  n  s  i  e  nt  cooling  produces 
stresses  in  a  grain  that  increase  with  time  (and  decreasing  temperature) 
toward  a  maximum  value.  The  maximum  stresses  are  attained  near  the  time 
that  thermal  equilibrium  of  the  grain  is  reached  and  decay  monotonically  as 
time  progresses.  The  strength  of  the  propellant  also  increases  as  the  tem¬ 
perature  decreases  and  will  be,  for  the  stated  conditions,  a  maximum  near 
the  thermal  equilibrium  time  of  the  motor.  At  longer  times,  the  comparative 
strength  will  decrease.  The  probability  of  failure,  as  evaluated  in  terms  of 
stress,  depends  upon  the  relative  rates  of  increase  of  grain  stress  and  grain 
strength  with  decreasing  temperature.  The  stress  rates  can  be  calculated  by 
the  methods  previously  discussed  (Section  4,  Vol.  I).  However,  the  strength 
increase  rate  with  temperature  cannot  be  interpreted  accurately  by  presently 
available  techniques  for  the  condition  of  continuously  varying  temperature. 

This  problem  area  merits  continued  study.  It  is  reasonable  to  expect  that 
grain  unbonding  as  a  result  of  thermal  stress  will  compare  to  a  propellant 
stress  criterion.  The  rationale  for  this  instance  arises  from  the  observation 
that  the  stresses  at  the  case-bond  interface  are  not  accompanied  by  signific¬ 
ant  straining,  except  near  the  ends  of  a  finite  grain. 

A  limit  condition  considered  in  interpretation  of  the  shock-heated  analog  motor 
failures  is  of  interest  (Section  5.  Z.  Z  above).  Viscoelastic  stresses  resulting 
from  imposition  of  the  thermal  shock  on  the  case  were  calculated  assuming 
no  change  in  the  grain  temperature.  The  increment  of  hoop  stress  at  the  port, 
introduced  by  the  case  expansion  exceeded  the  isothermal  stress  capability  of 
the  propellant  as  measured  in  the  laboratory  at  the  corresponding  strain  rate 
and  temperature.  The  case-bond  interface  stress  increment,  as  compared  to 
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the  propellant  strength  at  the  case  temperature,  also  exceeded  the  measured 
propellant  strength. 

However,  as  previously  noted,  the  incremental  port  hoop  strain  increment 
also  exceeded  the  propellant  strain  capacity  and  a  separation  of  criteria  was 
possible  only  in  that  port  cracking,  not  bond  failure,  was  observed. 
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APPENDIX  A 

THERMAL  FIELD  SOLUTIONS 


A-1.  INTERNAL  HEATING  OR  COOLING* 


A  hollow  cylinder  whose  outer  surface  is  insulated  and  whose  inner  port  is 
subjected  to  a  step  change  in  temperature  is  considered.  The  solution  as¬ 
sumes  a  temperature  distribution  as  follows: 


T  = 


(r-a) 


r-a  <  qj 
r-a  >  qj 


(A.l) 


As  shown  in  Figure  A- 1,  the  penetration  depth  qj  varies  with  time.  The  heat 
front  in  this  case  moves  from  the  constant  temperature  inner  port  to  the 
insulated  outer  surface.  Equation  A.  1  is  compatible  with  this  phenomenon 
and  is  therefore  a  suitable  choice  for  the  approximate  techniques  given  by 
Boley  and  Weiner,  I960,  The  relation  between  qj  and  the  time  was  obtained 
from  the  energy  technique  and  is  as  follows: 


^  -  30  Kt 

a*  ^  a'  '  a^ 

or 

^  +  3 qj  =  30X^  n" 

where 

qj  =  qj/a 

X  =  b/a 

=  Kt/b^ 


(A.  2) 


(A.  3) 


(A.  4) 


*  Reference  Section  2,  Volume  I  for  terminology. 
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Figure  A-1  Penetration  Depth  versus  Time 
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Again,  the  equation  for  the  penetration  depth  is  a  cubic  with  only  one  signifi¬ 
cant  root.  This  root  is  obtained  as  follows: 


q,  =  y-1 

and 

y  =  2  Cos  0/3 

where 

Cos  0  =  -  1  +  15T1^ 
Equation  A.  7  is  valid  only  for 


(A.  5) 


(A.  6) 


(A.  7) 


(A.  8) 


since 


Cos  0  <  1 


For  values  of  X<  2,  Equations  A.  5  through  A.  7  can  be  used  to  solve  for  qj 
for  each  value  of  and  hence  time.  However,  if  X.  5^2,  the  solution  of  the 
cubic  cannot  be  solved  in  this  simple  manner.  Since  qj  is  the  nondimensional 
penetration  depth,  physically  it  must  lie  between  1  and  X  -1.  Hence  the  value 
of  can  be  computed  for  each  value  of  time,  for  example,  by  an  itera¬ 
tion  technique. 

As  time  progresses,  q^  increases  until  it  obtains  the  value 


^13 


^  =  X-l 
a 


(A.  '1) 


or  the  heat  front  reaches  the  outer  surface.  The  time  interval  necessary  for 
this  to  occur  is  found  by  setting 

Tj  =  ^  1 

hence  from  Equation  (A.  3)  the  value  of  the  thermal  diffusion  parameter,  7]^ 
is  given  by 


-  3X 


+  3 


_  (X-l)"(X+2) 
- JOF - 


(A.  10) 
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Hence  the  time,  t^,  can  be  found  from  the  definition  of  (Equation  A.  4), 

It  is  interesting  to  note  that  as  ^ — >  1,  n^  - >  0;  and  as  X — >  * ,  — >  <®  . 

That  is,  the  heat  front  traverses  an  infinitely  thin  cylinder,  (X  =  l),  instan¬ 
taneously;  while  for  a  cylinder  with  a  "pin  hole"  for  an  inner  port,  (X  =“),  it 
takes  an  infinite  time.  It  is  also  interesting  to  note  that  the  expression  for 
(Equation  A- 10)  for  external  cooling  does  approach  a  finite  limit  as 
\ — >  ®.  That  is,  the  equations  for  a  hollow  cylinder  given  in  the  main  text 
are  also  valid  for  a  solid  cylinder.  This  is  accomplished  by  letting  X — > 

Once  the  temperature  front  reaches  the  outer  port,  a  new  function  must  be 
selected  in  lieu  of  Equation  (A.  1).  This  new  function  must  meet  the  following 
two  requirements: 

(1)  At  r  =  b,  the  temperature  must  equal  T^ 

(2)  Since  the  outer  surface  was  assumed  insulated,  the  derivative  of 
the  temperature  with  respect  to  "r"  must  equal  zero  at  r  =  b. 

Obviously,  Equation  (A.  1)  does  not  meet  these  requirements.  Hence  for  all 
t>  —  following  thermal  field  was  assumed: 

=  Ib:iF  I  (  ^'T^)  ^  T; 

Setting  r  =  b,  it  is  clear  that  T*  represents  the  temperature  at  the  outer 
insulated  surface  (Figure  A-2)  and  hence  is  a  function  of  time.  The  pro¬ 
cedure  for  finding  this  functional  relationship  is  similar  to  that  used  for  qj  . 
The  result  is  as  follows: 


=  T  [l-e''"* 


where 


(A.  12) 


(A.  13) 


_  (\  +  3)(\  +  2) 
1  =  HI  a  +  5) 


(A.  14) 
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An  inspection  of  Equation  (A.  12)  reveals  that  when 

or 

=  1 

then 

T4  =  0 


and  when 


Tl^  >«  ;  T4  >  Tq 

since,  as  stated  above,  T4  represents  the  temperature  at  the  outer  surface, 
the  relation  for  T4  is  physically  reasonable.  Hence,  an  approximate  solution 
for  this  case  results  from  the  use  of  Equations  (A.  1),  (A.  10)  and  (A.  11)  once 
qj  and  T4are  determined  from  Equations  (A.  3)  and  (A.  12).  The  results  for 
a  cylinder  with  a  b/ a  ratio  equal  to  2  are  shown  in  Figure  A-3, 


A- 2  INTERNAL  AND  EXTERNAL  HEATING  OR  COOLING 


A  hollow  cylinder  whose  outer  and  inner  surfaces  are  subjected  to  a  simul¬ 
taneous  step  temperature  change  is  considered.  This  solution  is  more  com¬ 
plicated  than  the  previous  two  because  instead  of  one  moving  heat  front  there 
are  now  two.  Hence,  it  is  necessary  to  choose  two  functions  for  the  tem¬ 
perature. 

The  assumed  functions  are: 


T,  (r.t)  .  { 

To  tl- 

r  -  a  <qj 

(A.  15) 

1 

^  0 

;  r  -  a  >  qj 

^2  (r,  t)  =  1 

To 

;  b  -  r  <  qj 

q. 

(A.  16) 

L  0 

;  b  -  r  S:  q 

A-6 


LOCKHEED  PROPULSION  COMPANY 


578-F 


A-7 


LOCKHEED  PROPULSION  COMPAN' 


Figure  A-3  Time -Temperature  Variation  in  a  Hollow  Cylinder  with  Zero  Initial 
Temperature  and  a  Suddenly  Applied  Surface  Temperature,  T^,  on 
the  Inner  Port  (outer  surface  insulated) 
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where  (Figure  A-4) 
T,  (r.t)  = 
T2(r.t)  = 


temperature  distribution  initiating  at  the  inner  port 

temperature  distribution  initiating  at  the  outer  surface 

penetration  depth  for  the  heat  front  proceeding  from  the 
inner  surface 

penetration  depth  for  the  front  moving  away  from  the 
outer  surface 


Once  again,  the  heat  fronts  progress  through  the  web.  Before  the  two  fronts 
meet,  material  points  lying  in  the  region  r  -  a  <  are  not  "aware"  of  the 
conditions  on  the  outer  surface.  Similarly,  material  points  in  the  region 
b  -  r  <  q,,  are  not  "aware"  of  the  conditions  on  the  inner  port.  Thus,  the 
variations  of  qj  and  q^  with  respect  to  time  are  the  same  as  those  given  in 
Section  A.  1  of  this  appendix  and  the  solution  in  the  main  text. 

Such  is  not  the  case  when  the  two  fronts  meet.  When  this  occurs,  the  two 
aforementioned  solutions  do  not  apply.  That  is,  one  cannot  superimpose 
those  temperature  distributions  to  obtain  the  distribution  for  the  present  case. 
Since  the  boundary  conditions  for  those  cases  are  not  equivalent  conditions 
for  this  case,  superposition  is  not  valid.  Thus,  the  values  of  the  respective 
penetration  depths  and  the  time  must  be  determined  at  the  instant  the  two 
fronts  meet.  This  is  accomplished  by  using  Equation  (2-10)  of  Section  2, 
Volume  I,  and  Equation  (A, 3)  of  the  present  section  as  follows: 


q’  +  3q^  -  30\^  Tl^  =  0 
q,  -  3q^  +  30T1"  =  0, 

where 

q3  =  ^3/  ® 

q,  =  qi/b  ■ 

and  from  the  geometrical  fact  that  when  the  fronts  meet 
q,  +  q3  =  b-a 
or 

q,  .  ^  (A.  19) 


(A.  17) 
(A.  18) 
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Figure  A -4  Time -Temperature  Distribution  of  a  Hollow  Cylinder  with 
Zero  Initial  Temperature  and  a  Suddenly  Applied  Surface 
Temperature,  T_,  on  the  Outer  and  Inner  Surfaces 
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The  simultaneous  solution  of  Equations  (A.  17,  18,  19)  for  the  unknown  quanti¬ 
ties  qjj^,  qjg^,  and  results  in  the  following: 

+  3  qja  (X+  1)  -  (2X+  1)  (X-1)  =  0  (A.  20) 


or 


or 


and 


-(^-1)  -  q,. 

^la 

K 

Tli 

=  "^a  +  3  q^ 

'a 

30X^ 

^.a* 

^sa*  ^a 

meet. 

For  example. 

^Ja 

=  0.47 

q 

=  0.235 

^3a 

^la 

=  0.265 

'lia 

=  0.265 

=  0.00639 

(A.  21) 
(A.  22) 


Thus,  the  fronts  meet  slightly  to  the  left  of  the  web  midpoint  (r  =  0.735). 


Once  the  fronts  meet,  two  new  functions  are  necessary  to  describe  the  tem¬ 
perature.  They  are  as  follows; 


li 

Tz 


(l-^  )  [r^-r]^  +  l^  [r^-a]^  | 

(1-^)  [ra-r]^  +  ^  [r^-b]^  } 


r  <  rg^  (A.  23) 

r  ^  (A.  24) 
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where 

=  distance  from  origin  at  which  the  two  fronts  meet 
T5(t)  =  temperature  at  r  =  r^  . 

The  dependence  of  Tj  with  respect  to  time  was  found  by  using  the  same  tech¬ 
nique  as  before.  The  expression  for  Tj  for  X  =  2  is 

Tj  =  Tq  [  1-e  ]  (A.  25) 

where 


Oj  =  0.2539 

The  thermal  distribution  for  a  cylinder  with  X  =  2  is  given,  for  various  times, 
in  Figure  A-4. 
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Appendix  B 

ON  THE  THEORY  OF  PHOTOVISCOELASTICITY* 


Abstract 


The  general  problem  of  photoviscoelasticity  cannot  be  solved  by 
a  simple  extension  of  the  experimental  and  interpretive  methods 
which  are  used  in  classical  photoelasticity.  For  viscoelastic 
materials  a  complete  time  history  of  isochromatics  and  isoclinics 
must  be  recorded.  The  stress  state  is  then  determined  by  an 
integral  operation. 


*  Ellis  Harold  Dill 

Department  of  Aeronautical  Engineering 
University  of  Washington 
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1 .  Summary 

The  stress  analysis  of  bodies  with  irregular  boundaries  is  such  a  complicated 
problem  that  it  is  frequently  desirable  to  determine  the  stress  state  by 
experiments  upon  a  model.  Photoelasticity  has  been  of  great  value  in  the 
design  with  elastic  materials,  and  it  would  be  desirable  to  have  this  tool 
available  in  all  cases.  However,  for  some  materials,  including  solid  propel¬ 
lants,  the  material  from  which  the  prototype  is  constructed  cannot  be  regard¬ 
ed  as  an  elastic  substance,  but  instead  must  be  treated  as  a  viscoelastic 
material.  In  general,  the  strain  state  in  the  viscoelastic  material  is  differ¬ 
ent  from  that  in  an  elastic  material  which  is  subjected  to  the  same  boundary 
conditions.  In  such  problems  the  viscoelastic  strains  cannot  be  obtained 
directly  from  observations  on  a  model  made  of  an  elastic  material.  This 
problem  of  correlation  between  model  and  prototype  would  be  merely  a  mat¬ 
ter  of  geometric  scaling  if  the  model  material  has  the  same  mechanical 
properties  as  the  prototype.  It  may  be  possible,  for  example,  to  manufacture 
a  model  of  the  desired  viscoelastic  properties  from  a  high  polymer  by  proper 
selection  of  epoxy  resin,  plasticizer,  and  hardener. 

Photo-analysis  is  based  upon  the  observation  that  certain  materials  exhibit 
polarization  which  depends  upon  the  mechanical  state  variables.  Such  materi¬ 
als  are  termed  birefringent.  They  are  electromagnetically  anisotropic, 
except  perhaps  in  the  natural  state.  The  propagation  of  a  plane  electro¬ 
magnetic  wave  through  a  slice  of  such  material  is  reviewed  in  Section  2.  It 
is  shown  that  the  wave  is  completely  linearly  polarized  in  certain  planes. 

Thus,  an  incident  wave  is  divided  into  two  waves  propagating  at  different 
velocities  characterized  by  the  principal  indices  of  refraction.  By  observa¬ 
tion  on  the  emanating  wave  it  is  possible  to  determine  the  refractive  indices 
and  thus  the  dielectric  constants. 

The  principles  of  linear  thermoviscoelasticity  are  reviewed  in  Section  3.  For 
linear  viscoelastic  materials,  stress  and  strain  are  uniquely  related  to  each 
other  by  integral  or  differential  operators. 

A  viscoelastic  material  may  be  expected  to  exhibit  optical  as  well  as  mech¬ 
anical  time  effects,  such  as  creep  at  constant  stress.  The  dielectric  con¬ 
stants  which  determine  the  refractive  indices  will  depend  not  only  on  the 
instantaneous  values  of  the  state  variables,  but  also  on  the  history  of  mech¬ 
anical  states.  In  the  case  of  linearly  viscoelastic  materials,  the  mechanical 
state  variables  may  be  either  stress  or  strain.  The  one  being  related  to  the 
other  uniquely,  they  are  interchangeable.  For  the  small  changes  in  the 
mechanical  state  variables  which  occur  in  linearly  viscoelastic  materials,  it 
would  be  expected  that  a  linear  relation  between  the  dielectric  properties  and 
the  mechanical  state  variable  could  be  used.  The  determination  of  the  stress 
state  from  observations  on  the  birefringence  of  materials  for  which  such  a 
relation  exist  is  discussed  in  Section  4.  It  is  found  that  the  principle  axes  of 
stress  do  not,  in  general,  coincide  with  the  principal  axis  of  polarization  as 
they  do  in  elastic  materials.  Therefore,  the  difference  in  principal  stresses 


B-2 


LOCKHEED  PROPULSION  COMPANY 


578-F 


is  not  determined  solely  by  the  difference  in  the  indices  of  refraction,  but 
also  by  the  angle  between  the  principal  axes  of  polarization  and  the  principal 
axes  of  stress.  Thus,  for  practical  stress  analysis,  routine  methods  of 
recording  both  the  difference  in  indices  of  refraction  (relative  retardation) 
and  the  plane  of  polarization  (isoclinics)  as  a  function  of  time  must  be  devised. 

In  the  following  sections,  all  quantities  are  referred  to  a  stationary  rectangu¬ 
lar  cartesian  coordinate  system  (z*,  z^,  z*)  =  (x,  y,  z).  Vectors  are  denoted 
by  subposed  tilde.  The  time  is  denoted  by  t.  The  symbols  are  defined  as 
they  are  encountered.  Brackets  are  used  in  Section  2  to  indicate  the  difference 
in  a  quantity  on  two  sides  of  a  surface.  The  base  vectors  are  denoted  by  g.. 

All  indices  have  the  range  1,  2,  3.  Repeated  indices  indicate  summation.** 

The  symbol  [  g,,  fe,  g]  indicates  the  scalar  triple  product  of  the  vectors.  In 
Section  2,  the  dielectric  tensor,  although  dependent  on  time,  is  assumed  to 
change  slowly  compared  to  the  wave  speeds  and  thus  can  be  regarded  as 
constant.  The  velocity  vector  is  y. 


2.  Electromagnetic  Waves 
2. 1  Fundamental  Equations* 

The  principles  of  conservation  of  charge  and  magnetic  flux  lead  to  the  follow 
ing  field  equations  which  hold,  except  on  certain  singular  surfaces: 


div  J  +  —  =  0,  (2.01) 

3t 

curl  E  +  i-i-  =  0,  (2.02) 

dt 

div  B  =  0.  (2.03) 

Here,  B  is  the  density  of  magnetic  flux,  E  is  the  electric  field,  J  is  the 
current  density,  and  Q  is  the  charge  density. 

These  field  equations  imply  the  existence  of  electromagnetic  potentials  such 
that 


B  =  curl  A,  (2.04) 

E  =  -1^  -  grad  V,  (2;05) 

dt 


*  References  are  listed  at  the  end  of  the  appendix 
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Q  =  div  p,  (2.06) 

J  =  curl  H  -  .  (2.07) 

Here,  A  is  the  magnetic  potential:  V  is  the  electric  potential;  H  is  the  cur¬ 
rent  potential  (usually  called  the  magnetic  field  intensity);  D  is~the  charge 
potential  (usually  called  the  electric  displacement). 

The  electric  field  and  the  magnetic  flux  density  are  related  to  the  charge 
potential  and  the  current  potential  at  all  points,  inside  and  outside  of  the 
material,  by  the  "aether"  relations: 

p  =  eo  E,  (2.08) 

H  r  J-  B.  (2.09) 

1^0  ~ 

Cg  and  p,g  are  fundamental  constants  depending  on  the  units  of  time,  charge, 
and  magnetic  flux. 

Across  a  surface  of  discontinuity,  such  as  the  boundary  between  regions  of 
different  material  properties,  which  moves  with  speed  Uj^  and  has  the  unit 
normal  v  ,  the  following  jump  conditions  hold: 


V  X 

[E  ]  -  u„  [b]  =  0, 

(2.10) 

[b]  •  V  =  0, 

(2.11) 

[J] 

•  -  ^n  tQ  ]  =  0. 

(2.12) 

In  addition  to  these  fundamental  fields  and  potentials,  constitutive  equations 
defining  the  electromagnetic  nature  of  the  material  are  required.  The  formu¬ 
lation  of  constitutive  equations  is  facilitated  by  the  introduction  of  certain 
auxiliary  fields  and  potentials:  I.,et  P  be  the  density  of  polarization  and  M 
the  magnetization  field.  The  polarization  current  Jp,  the  bound  charge 
Qg,  and  the  magnetization  current  J  are  given  by 


ip 

dp 

=  +  curl  (P  X  V  ), 

(2.13) 

Qb 

"  -  fliv  P. 

(2.14) 

iM 

=  curl  M  . 

(2.15) 
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The  total  charge  is  considered  to  be  made  up  of  the  bound  charge  plus  the 
free  charge: 

Q  =  Qfi  +  Qf 

The  total  current  is  the  sum  of  the  free,  the  polarization,  and  the  magneti¬ 
zation  currents: 

i  ^  +  ip  +  iM-  (2.17) 

In  terms  of  the  partial  potentials  H  and  D  which  are  defined  by 

=  H  -  M  -  P  X  V  ,  (2,18) 

p'  =  p  +  P  ;  (2.19) 

the  free  current  and  free  charge  are  given  by 

Jp  =  curl  h'  -  ^  ,  (2.20) 

~  a  t 

Qp  =  div  d'  (2.21) 

The  jump  conditions  across  a  moving  surface  of  discontinuity,  with  speed  un 
and  normal  v  are 

[p'l  •  V  =  0,  (2.22) 

V  X  Ch]  +  Ujj  [p5  =  0.  (2.23) 

2, 2  Ideal  Dielectrics 

An  ideal  anisotropic  dielectric  is  defined  by  the  equations 
J  F  =  0, 

Qf  -  0,  (2.24) 

M  =  0, 
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Therefore 


where 


Di 


=  e 


0  Cj, 


K 


ij 


=  6 


ij 


+  Pii 


(2.25) 


(2.26) 


The  fundamental  equations  for  dielectrics  are  equations  (2.25),  (2.02),  (2.03), 
and  (2.09),  combined  with  equations  (2,20)  and  (2,21)  which  become 


curl  H 


d  d'_ 


0, 


(2.27) 


div  D  =  0. 


(2.28) 


The  K-  will  be  called  the  dielectric  tensor. 


2.  3  Plane  Electromagnetic  Wave 


Suppose  that  a  dielectric  at  rest  is  transversed  by  a  plane  electromatnetic 
wave  given  by 


E  =  e 

PM  M 

ei'P, 

(2.29) 

B  =  b 

ei'P, 

(2.30) 

p'  =  d 

ei'P. 

(2.31) 

h'  =  h 

ei'P, 

(2.32) 

9  =  k 

n  •  z  -  uJt. 

(2.33) 

The  unit  vector  n  gives  the  direction  of  the  wave;  k  is  the  wave  number;  is 
the  angular  frequency;  e,  b,  d,  and  hare  constant,  complex  vectors. 


Substituting  equations  (2.29)  -  (2.33)  into  equations  (2.02)  and  (2.27)  gives 


knxe=bu),  (2. 34) 

k  n  X  h  +  d  =  0.  (2.35) 


Therefore,  b  is  perpendicular  to  n  and  e  while  d  is  perpendicular  to 
Equation  (2.59)  shows  that  b  is  parallel  to  h  for” a  dielectric  at  rest, 
n,  d,  h  are  mutually  perpendicular  and  d,  e,  n  lie  in  the  same  plane 
ure"fe-T).  ~  ~  ~ 


n  and 

~Thus 

(Fig- 
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Figure  B-1  General  Plane  Wave 


The  plane  d,  h  is  the  wave  front.  The  direction  of  h  is  the  direction  of 
polarization.  The  plane  n,  h  is  the  plane  of  polarization.  The  direction  of 
is  the  direction  of  propagation.  The  direction  of  e  is  the  direction  of  vibra( 


The  constitutive  relation  for  a  dielectric,  equation  (2.25),  gives 

.1 

^  ®i  ”  ^ij 


(2.36) 


Equation  (2,36)  gives 


-1 


(uSq  ®jrs  ^s’ 


(2.37) 


Equations  (2.09)  and  (2.35)  lead  to 


-1 


(6ik  +  N  e.pq  n^)  hj^  =  0. 


(2.38) 


Where, 


is  the  index  of  refraction  and 


c  = 


is  the  speed  of  light. 


(2.39) 


(2.40) 
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Let  us  choose  to  be  the  principal  axes  of  Kjj  and  consider  the  case  when  n 
is  directed  along  z.  Then  equation  (2.38)  becomes  ~ 


1  -  nVKz 
0 
0 


1  -  nVk, 
0 


1 


=  0 


(2.41) 

Equations  (2.41)  and  2.37)  show  that  only  two  plane  waves  are  possible;  Either 


or 


N  =  Ni  = 


hi 


h3 


62  =  ej 


0,  hz  =  N]  Sp  ei, 


(2.42) 


N  =  Nz 


hz  =  hs  =  ei  =  ej  =  0,  hi  =  -Nz  ez. 


(2.43) 


The  direction  of  polarization  coincides  with  one  of  the  principal  axes  and  the 
direction  of  vibration  is  normal  to  the  wave  direction  (Figure  B-2).  Every 
wave  must  be  a  linear  combination  of  these  two  waves. 


FigureB-2  Plane  Wave  Propagating  Along  One  Principal  Axis,  Case  1 


In  free  space  =  6  —  ;  therefore  N  =  1  and  hs  =  63  =  0. 


2.  4  Plane  Polariscope 

Suppose  a  slice  of  dielectric  is  bounded  by  the  planes  z  =  0,  h.  Consider  the 
case  when  the  coordinate  axes  are  principal  axes  of  the  dielectric  tensor  Ky. 
By  means  of  a  polarizer,  a  linearly  polarized  plane  electromagnetic  wave 
propagating  along  the  z  axis  is  obtained.  Every  such  wave  can  be  divided  into 
the  sum  of  monochromatic  waves  and  every  monochromatic  wave  can  be 
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decomposed  into  components,  given  by  equations  (2,29)  -  (2.33),  such  that 
the  direction  of  polarizarion  of  each  component  coincides  with  one  of  the 
coordinate  axes.  Let  this  plane  wave  be  incident  on  the  surface  z  =  0.  In 
order  to  satisfy  the  jump  conditions,  equations  (2.10)  and  (2.23),  there  must 
also  be  a  wave  traveling  in  the  opposite  direction  (reflected). 


The  electromagnetic  field  preceding  the  dielectric,  Figure  B-3,  is  therefore 
the  sum  of  two  fields: 


E  = 


~(1)  ^  ~  (2) 


(2.44) 


where 


z-tut) 

E(l)® 


z-U)t) 


The  wave  transmitted  into  the  material  is  itself  reflected  at  the  surface  z=h. 
The  electromagnetic  field  in  the  dielectric  material  is  therefore 

?  ^  ~(3)  ~  (4)  (2.45) 


where 


E 


(4)  =  S(4) 


i(-kz 

e 


-UJt) 


The  wave  transmitted  across  the  surface  z=h  gives  rise  to  the  electromag¬ 
netic  field  following  the  dielectric: 


where 


~(5) 


^i(ko  z-u)t) 


ho  - 


U) 

c 


(2.46) 
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N(5) 


1 

E 


FREE  SPACE 


Figure  B-3  Directions  of  Plane  Waves 


The  vectors  e/.  v  are  complex  constant  vectors  which,  in  the  present  case,  are 
all  directed  along  the  same  one  of  the  coordinate  lines.  The  value  of  k  used 
in  equation  (2.45)  is  determined  from  the  appropriate  one  of  equations  (2.42) 
or  (2.43),  The  corresponding  magnetic  field  is  determined  from  equation 
(2.34). 

The  jump  conditions,  applied  to  the  plane  z  =  0  and  to  the  plane  z  =  h,  lead  to 
a  set  of  four  vector  equations  containing  the  five  e,, ,  which  yield  the  following 
result*  ~  ' 


^(5) 


4N  e  ** 


(1+N)*  -(1-N))*  e* 

If  N  is  of  order  of  magnitude  unity,  then  approximately* 

.  4N 


(1+N)* 

The  further  simplification* 


i(k-ko)h  ■  ’ 

®  E(l)* 


(2.47) 


(2.48) 


•  i(k-ko)  h 

£(5)  =  ®  £(i) 


is  also  sometimes  a.isumed. 


(2.49) 


By  assigning  the  appropriate  value  to  k,  this  result  applies  to  incident  waves 
polarized  along  either  of  the  coordinate  axes. 
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Consider  an  incident  monochromatic  transverse  wave  obtained  by  passing 
monochromatic  light  through  a  polarizer: 

Ej  =  (A  cos  O'  g,  +  A  sin  » gz )  (2.50) 

where  a  is  the  angle  between  Ej  and  the  x  -  axis. 

Applying  equation  (2,47)  to  each  component  of  the  incident  wave,  the  trans¬ 
mitted  wave  is  found  to  be 


E-p  =  A(ci  cos  ogj  +  cj  sin  ot  g^)  (2.51) 


where 


C|  = 


-ik„h 
4Ni  e  0 


(1+N,)' :(1-N,)^ 


^  _  4N2e'^‘^o^ 

(l+Nz)^  -(l-Nz)^  e‘^^^ 

Passing  this  wave  through  a  second  polarizer,  will  give  a  wave  of  complex 
amplitude 

E  =  A(c,  cos  a  cos  3  +  c^  sin  a  sin  3  )  (2.52) 

where  3  is  the  angle  between  the  transmitted  wave  and  the  x  -  axis. 

This  can  be  put  in  the  form 

E  .  4  A  K:  cos  «  cos  3  +  9,)  ^^  33^ 

,  A  A  V  ■  n,  ■  R  i(koZ-<"t-koh  +  9  z) 

+  4  A  Kz  sin  O'  sin  P  e  ” 

where 

Kn  = -  (2.54) 

n  -  z 

2v  4  cos  kj,h  +  (l  +  N^)  sin^  k^h 
tan  tan  k^^h.  (2.55) 
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Either  the  real  or  imaginary  part  of  E,  equation  (2.53),  may  be  regarded  as 
the  light  transmitted.  The  imaginary  part  is 

I  =  Q  sin  (koz  -  tut  -  koh  -  x)  (2.56) 


where 


Q  cos  X  =  4  A  K|  cos  o  cos  P  sin  6  j 

+  4  A  K2  sin  ot  sin  P  sin  6  j  , 

Q  sin  X  =  4  A  Ki  cos  a  cos  P  cos  6  1 

+  4  A  Kz  sin  a  sin  P  cos  P  2  . 

Consider  the  case  P  =  a  +  n / 2,  i.  e. ,  the  first  and  second  polarizers  are 
crossed.  Then 


=  4  A^  sin^  a  [  K,^  +  Kz^  -  2Ki  Kz  cos  X  n  „  ]  (2.57) 


where  the  fringe  order  n  is  defined  by 

(2.58) 

Therefore  the  light  intensity  is  zero  where  a  =  0,  n/ 2.  Lines  connecting  all 
such  points  are  called  isoclinic  lines.  The  light  intensity  will  be  a  local 
minimum  where  Ni/Nz  =  Cj,  Cz.  The  lines  connecting  all  such  points 

are  called  isochromatics. 


n  = 


62-9 

2n 


For  the  transparent  materials  used  in  photoanalysis  Ni  is  equal  to  Nz  in  the 
natural  state,  and  both  Nj  and  Nz  change  very  little  as  a  result  of  deformation 
of  the  material.  In  such  a  case,  the  minimum  values  of  equation  (2.57)  occur 
at  integer  values  of  the  fringe  order  n.  If  Nj  and  Nz  are  sufficiently  close  to 
unity,  then  equation  (2.55)  gives 


n 


(Nz  -  N,). 


(2.59) 


Therefore,  for  materials  with  an  index  of  refraction  which  is  near  unity  and 
which  changes  only  slightly  upon  deformation,  the  isochromatic  lines  occur 
at  definite  values  of  Nz  -  Ni  as  given  by  equation  (2.59). 
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3.  Viscoeiatiticity 

3.  1  Fundamental  Equations 

The  principle  of  conservation  of  mass  leads  to  the  field  equation 

I-?  +  div  (pv)  =  0  (3.01) 

o  t  ~ 

where  p  is  the  density  and  v  is  the  velocity  vector. 

The  stress  vector  t  on  a  surface  with  unit  normal  n  is  related  to  the  stress 
tensor  t^ni  by  the  Equation  ~ 

tk,n)=  (3.02) 

The  principle  of  conservation  of  linear  and  angular  momentum  leads  to  the 
equations  of  motion: 

m  +  P  fl'  =  p  a^  (3.03) 

tkm  =  (3.(H) 

where  f  is  the  body  force  and  a  is  the  acceleration. 

In  generalf  the  electromagnetic  field  gives  rise  to  the  bory  force 

f  =  QE  +  JxB. 

«*«««#«« 


3.  2  Constitutive  Equations  for  Deformation^’  ^ 
The  deviatoric  stress  a^j  is  defined  by 

"  *ij  ■  T  ‘kk  *  ij- 

The  quantity 

P  =  -  T  ‘kk 
is  the  mean  pressure. 


(3.05) 


(3.06) 
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The  deviatoric  strain  is  defined  by 

®  '  T  ^ij*  (3.07) 

For  small  strains,  the  relative  volume  dilatation  is 


dv  -  dV 
dV 


®kk  =  e. 


(3.08) 


Empirical  evidence  suggests  that,  for  some  materials  experiencing  isotherm¬ 
al  deformation  and  small  displacements,  the  material  is  characterized  with 
sufficient  accuracy  by  the  relations: 


dx 


J(t-x)  dx 


(3.09) 


dx 


The  inverse  of  equations  (3.09)  and  (3.10)  has  the  form 


<7ii  =  J.„  G(t-x) 


d  ^^.(x) 

dx 

p  =  -  J.«,K(t-x)  ^5^  dx 


dx 


(3.10) 


(3.11) 

(3.12) 


Any  material  which  can  be  adequately  characterized  by  equations  (3.09, 

3. 10)  or  by  (3. 1  1,  3. 12)  is  called  "linear  viscoelastic.  "  The  kernels  of  these 
Duhamel  integrals  have  the  names 

G(t)  relaxation  modulus  in  shear. 

J(t)  creep  compliance  in  shear. 

K(t)  bulk  relaxation  modulus. 

M  (t)  bulk  creep  compliance. 

Alternative  notation  is 


G(t)  -  Go  +  ^  (t), 

J(t)  =  Jo  +  t/Tl  +  'l'(t), 
etc. 
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where  Jp  is  the  initial  elastic  compliance. 

"n  is  a  Newtonian  viscosity  coefficient. 

7  is  the  creep  function,  'f  (O)  =  0. 

Go  is  the  equilibrium  elastic  modulus. 

$  is  the  relaxation  function,  #  (“)  =  0 

This  notation  has  the  advantage  that  the  creep  function  is  generally  found  to 
be  a  monotonically  increasing  function  with  initial  value  zero  and  having 
monotonically  decreasing  slope.  The  relaxation  function  is  then  a  mono¬ 
tonically  decreasing  function  tending  to  zero  and  having  monotonically 
increasing  (algebraically)  slope. 

Each  of  these  quantities  is  temperature  dependent  as  well  as  time  dependent 
so  that  they  should  be  written  in  the  form  G  (t,  8),  where  6  is  the  absolute 
temperature.  For  isothermal  deformations,  empirical  evidence  suggests 
that,  for  many  materials,  the  temperature  dependence  can  be  adequately 
described  with  the  aid  of  a  "temperature  shift  function.  "  It  is  found  that, 
for  test  at  any  given  temperature,  if  a  reduced  time  T  defined  by 

InT  =  In  t  -  f(9)  (3.13) 

is  introduced;  then 

G(t(T))  =  G  (t),  (3.14) 

where  G  (t)  is  the  relaxation  modulus  at  a  reference  temperature  6  o-  The 
function  f(0  )  is  a  monotonically  decreasing  function  of  temperature  such 
that  f(9o)  =  0  and  is  regarded  as  a  material  property.  The  shift  factor  a(6) 
is  defined  by 


In  a  =  f(9  )- 


Then 


t  =  a(e)  T  (3.15) 

Similar  relations  are  assumed  to  hold  for  the  bulk  modulus.  A  material 
characterized  by  equation  (3.14)  is  said  to  be  thermo- rheologically  simple. 

Equation  (3.14)  implies  that 

G(0)  =  G(0)  =  Go  +  $(0) 

and  _ 

G(«')  =  G(»)  =  Go. 
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That  is,  the  initial  elastic  modulus  and  final  modulus  are  not  affected  by 
temperature  changes.  Thus,  only  the  relaxation  function  i  (t)  need  be  con¬ 
sidered: 


§(t)  =  $(T)  (3.16) 

where  $  (t)  is  the  relaxation  function  at  the  reference  temperature  Qq" 
Similarly,  for  the  creep  compliance: 

J(t)=J(T)  (3.17) 

implies  Jo  and  T)  are  independent  of  temperature  and 

'l'(t)  =  T(t)  (3.18) 

where  ^  (t)  is  the  creep  function  at  the  reference  temperature  6  q. 


3,  3  Constitutive  Equations  for  Linear  Thermoviscoelasticity 

There  is  a  viscous  dissipation  of  energy  but  if  the  heat  due  to  this  dissipation 
can  be  neglected  and  since  the  volume  changes  are  small,  the  temperature 
distribution  is  determined  by  the  uncoupled  Fourier  heat  conduction  equation. 

^  °  (3.19) 

Thus  the  temperature  may  be  regarded  as  a  known  function  of  space  and 
time  0=  0(x,  t).  The  theory  is  then  rendered  complete  by  specifying  a  phe¬ 
nomenological  relation  between  stress,  strain,  and  temperature. 

The  shear  creep  compliance  measured  by  a  creep  test  at  varying  tempera¬ 
ture  is  f(t).  If  the  initial  and  final  compliance  for  a  material  not  experiencing 
viscous  flow  are  independent  of  temperature,  then  f(0)  =T(0)and  f(“)  =  J(®), 
where  J(t)  is  the  creep  compliance  at  constant  reference  temperature  0o.  In 
this  case  there  is  a  function  §  (t)  called  the  reduced  time  such  that 

f(t)=J(Ut))  (3.20) 

and 

f'(t)  =  ^J'(§)  (3.21) 

where  the  prime  indicates  differentiation  with  respect  to  the  argument. 
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The  isothermal  creep  test  of  a  thermo -rheologically  simple  material  (equa¬ 
tion  3. 17)  yields 

j'(t)  =  a(0)  J'(T).  (3.22) 

It  is  physically  appealing  to  assume  that  this  relation^  holds  also  for  the 
creep  at  non-constant  temperature.  With  this  fundamental  hypothesis 

af  =  a(e(t)).  (3.23) 


Therefore 


(3.24) 


Now  suppose  that  the  same  temperature  field  exists,  but  the  load  is  applied 
at  time  t.  The  resulting  creep  compliance  is  F(t,  t  ):  The  total  strain 
resulting  from  a  continuous  distribution  of  load  is  then 

t  da.j(x) 

"  ifl  - ar -  (3.25) 


With  assumption  F(t,  t  )  =  j(0)  and  F(®,  t  )  =  J{“),  there  exists  a  function 
§  (t,  t)  such  that 


F(t,T)  =  J(§). 

(3.26) 

and 

dF(t,  T  )  _  d5  T-'/ir  \ 

—at - ar 

(3.26a) 

According  to  the  fundamental  hypothesis 

^  =  .<9(.)). 

(3.27) 

Since  F(t  , 

T )  =  J(0),  §  =  0  at  t  =  T . 

Therefore 

5(t.T)  =  a(9(Tl)  dT) 

(3.28) 

=  §(t)  -  l(T) 

and 

F(t.T)  =  J(5(t)  -  §(T)). 

(3.29) 

B-17 


LOCKHEED  PROPULSION  COMPANY 


578-F 


Introducing  the  notation 

eij(5(t))  =  eij(t)  (3.30) 

CTij(  (t))  =  Oij(t)  (3.31) 

the  constitutive  equation  becomes 

?  daij  (T) ) 

Sij  (?)  =  Jo  J(5-n)  — -  dq  (3.32) 

The  fundamental  constitutive  equation  for  thermoviscoelasticity  is  therefore 
equation  (3.25)  which,  under  the  fundamental  hypothesis,  can  be  put  into  the 
form  (3.32)  for  thermo-rheologically  simple  materials. 

4.  Photothermoviscoelasticity 

4.  1  Dielectric  Constitutive  Equations  for  Isothermal  Deformations 

Experimental  evidence  suggests  that  certain  polymers  exhibit  polarization 
under  isothermal  deformations  which  depend  upon  the  history  of  deformation. 
This  suggests  that  a  relationship  exists  between  the  dielectric  tensor  and.  the 
strain  tensor  similar  to  that  for  linear  viscoelasticity. 

Let  Nij  be  a  tensor  called  the  refraction  tensor  whose  principal  axes  coincide 


with  those  of  the  dielectric  tensor  Kjj,  and  such  that  the  principal  values 

Nj  of  N^j  are  related  to  the  principal  values  Kj  of  Kjj  by 

Ni  =VKi- 

(4.01) 

Denote  the  deviatoric  part  of  N^j  by  Njj  and  let 

N  .  i  Nkk- 

Consider  the  hypothesis 

^ij=J  X(t-T)^  [Sij(T)]dT, 

(4.02) 

N  =  J  A(t-T)  ^  Ce(T)]  dT. 

(4.03) 

where  x(t)  is  the  optical -relaxation  function, 

A  (t)  is  the  optical -volume -relaxation  function. 
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Denote  the  Laplace  transform  of  a  function  by  a  superposed  asterisk.  Then 


ITfj  =  p  X*  ejfj 

where  p  is  the  transform  variable.  From  equation  (3.13) 

e*  =  pj*  o  *  . 

ij  ^  ij 

Therefore 

a*  =  P  r*  N?., 
ij  ij 

where 

r*  =  — 1 - . 

pJ  X*  J”!” 

The  inverse  transform  gives 

a.j(t)  =  J„'r(t-T)^  [FIij(T)]dT. 


(4.04) 


(4.05) 


(4.06) 


(4.07) 


Thus,  a  unique  relation  between  the  polarization  and  the  history  of  deforma¬ 
tion  implies  a  unique  relation  between  stress  and  the  history  of  polarization 
and  inversely. 

Both  the  function  x(t)  and  F  (t)  may  be  expected  to  depend  on  the  temperature. 
Analogous  to  the  situation  for  the  mechanical  state  the  existence  of  an  optical- 
temperature  shift  factor  b(0)  may  be  supposed  such  that 

X(t)  =  X  (b  (0)t)  (4.08) 

where  x(t)  is  the  optical -relaxation  function  at  temperature  0q  and  x(t)  is  the 
optical-relaxation  function  at  temperature  0. 

Similarly 


r(t)  =  r  (b  (0)t). 


(4.09) 
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4.  2  Dielectric  Constitutive  Equations  for  Thermal  Stress 

Consider  deformations  that  occur  in  the  presence  of  varying  temperatures 
0(t).  If  the  polarization  at  various  strain  levels  but  the  same  temperature 
history  can  be  superimposed,  the  constitutive  equation  for  the  isothermal 
deformations  (4.02)  becomes 

N..  =  f(t.T)  ^  [  e..(T)]dT.  (4.10) 

The  function  f(t,  t  )  may  in  general  depend  upon  the  history  of  temperature. 

In  the  case  of  thermo- rheologically  simple  materials,  it  may  be  that  a 
reduced  time  p  (t)  exists,  in  general  different  from  that  for  mechanical  case, 
such  that 


f(t,  t)  =  X  (p(t,  t)). 

(4.11) 

t 

P(t,  t)  =  J_  b(0  (t))  dt. 

(4.12) 

In  this  event,  the  optical  properties  are  completely  characterized  by  optical 
relaxation  function  at  reference  temperature  9,  X(t),  and  the  optical- 
temperature  shift  factor  b(6).  Equation  (4.10)  can  then  be  written  in  a  form 
similar  to  equation  (3.36); 

Nij  (V  =1^  X  (T1-t)  [€..(t)]  dT.  (4.13) 


where 


T1(t)  =  J' 


t 


b(9  (t)  )  dt. 


Nij(Tl)  =  N.j  (t  (Tl)  ), 
(Tl)  =  ¥..(t  (Tl)  ). 


If  the  temperature  shift  factors  a(6)  and  b(0)  are  different,  no  simple  trans¬ 
formation  to  the  form  analogous  to  equation  (4.07)  exist.  But  if  a(0)  =  b(0), 
then  Tl  (t)  =  5  (t).  Taking  the  Laplace  transform  with  respect  to  5  of  equation 
(4.13)  and  equation  (3.36)  gives 


€. . 

ij 


-  * 

px 


(4.14) 

(4.15) 
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Therefore 


where 


5*.  =  p  r*N?., 


The  inverse  transformation  gives 


r(?-T)^  [N..  (T)]dT. 


(4.16) 


(4.17) 


4.  3  Principal  Axes 

Let  be  the  fixed  rectangular  cartesian  axis  with  -  z.  Consider  the  plane 
stress  problem  and  take  ^  =  z  normal  to  the  plane  of  the  .stress  field.  Sup¬ 
pose  that  z  is  a  principal  axis  of  a..,  e  . and  N^..  Let  x^  be  the  principal 
axis  of  strain,  x^  =  z*,  and  let  Of  be*  the^ingle  from  z*  to  x*.  Let  be  the 
principal  axis  of  stress,  =  z*,  and  let  3  be  the  angle  from  z*  to  y*.  Let 

x^  be  the  principal  axis  of  N^j,  x^  =  z^,  and  y  be  the  angle  from  z*  to  x*. 


Figure  B-1  Principal  Axes 


In  a  manner  similar  to  that  of  Section  3.  4,  it  follows  from  equation  (4.07)  that 


Ou(t)  -0  22(t)  -  r(t-T)^  [(N,{t) 

^I2(t)  =  r  (t-T  )  ^  C  (N,(T  )  -  N2(T  ) 


N2(t  )  )  cos 
sin  2Y(T  )] 


2Y(t)]  dT, 

(4.18) 


dT  . 


(4.19) 
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Using  equations  (3.32)  and  (3.33)^  it  is  seen  that  a  knowledge  of  the  angle 
Y  (t)  and  the  difference  Ni(t)  -  NzCt),  of  the  principal  values  of  N. as  a 
function  of  time  is  necessary  and  sufficient  to  determine  the  difJirence  in 
principal  stresses  and  the  principal  axes  of  stress  in  the  isothermal  case. 

A  similar  reasoning  applies  to  the  transient  temperature  case  when  equation 
(4. 17)  may  be  used. 


4.  4  Experimental  Procedure 

The  model  is  placed  in  a  standard  plane  polariscope  and  subjected  to  a 
desired  loading  and  heating  sequence.  The  isoclinic  and  isochromatic  patterns 
will  vary  with  time.  A  complete  record  of  the  time  variation  of  the  fringe 
order  n(t)  and  the  axis  of  polarization  y  (t)  at  a  point  is  made.  This  will 
require  a  more  elaborate  recording  device  than  is  normally  used  with  the 
polariscope.  For  a  given  model  thickness  h,  equations  (2.59)  and  (2.55) 
determine  a  value  for  Ni(t)  -  N2(t).  Equations  (4.18)  and  (4.19)  may  then  be 
used  in  the  isothermal  case  to  determine  the  difference  in  principal  stresses 
<7i(t)  -  Oi{t)  and  the  axis  of  principal  stress  3  (t)  as  functions  of  time.  To 
determine  the  material  optical  properties,  this  procedure  is  applied  to  any 
configuration  for  which  the  stress  is  known;  for  example,  a  creep  test. 
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